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1.  INTRODUCTION 


The  discovery  of  optical  second  harmonic  generation  (SHG)  in  1961,  coupled  with  more  recent 
advances  in  laser  technology  and  nonlinear  materials,  has  led  to  increasingly  rapid  progress  in 
the  application  of  nonlinear  materials  to  nonlinear  frequency  conversion  [1,  2,  3,  4]. 

This  document  reports  the  results  in  the  investigation  of  three-wave  interactions  in  positive 
biaxial  (nonlinear)  crystals.  The  results  are  characterized  in  terms  of: 

•  Phase  matching  angles, 

•  Effective  nonlinear  coefficients,  and 

•  Walkoff  angles  for  the  optical  beams. 

Specifically,  this  report  concentrates  on  KTi0P04  (KTP)  and  its  crystallographic  isomorphs: 
RbTi0P04  (RTP),  Kl'i0As04  (KTA),  RbTi0As04  (RTA),  and  CsTi0As04  (CTA).  The 
isomorphs  KTA,  RTA,  CTA,  and  RTP  are  relatively  new  and  show  promise  as  nonlinear  crystals. 
Potassium  Titanyl  Phosphate  (KTP)  is  a  widely  used  nonlinear  optical  material  which,  since  there 
is  a  large  library  of  information  for  this  particular  crystal,  was  used  to  validate  the  computer 
model  calculations  developed  for  this  study. 

Prior  to  this  study,  most  calculations  assumed  that  the  electric  field  S  was  perpendicular  to 
the  wave  vector  I?  and,  therefore,  replaced  ^  by  the  electric  displacement  0.  Also,  when  the 
walkuff  angle  was  calculated,  the  fact  that  optical  waves  polarized  in  different  directions  have 
different  walkoff  angles  was  ignored.  The  consequence  was  that  the  approximate  treatment  was 
prone  to  the  introduction  of  errors  or  inaccuracies  in  the  final  results. 

This  report  will  show  the  results  of  a  translatioti  of  an  exact  mathematical  theory  developed 
by  Yao  and  Sheng  [3]  into  a  practical  computer-based,  mathematical  model.  The  model  is  of  a 
generalized  nature  which  uses  accurate,  rather  than  approximate,  calculations. 
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2.  ANGULAR  PHASE  MATCHING 


a.  Phase  Matching  Theory. 

Nonlinear  optica!  interactions  in  bulk  materials  (e.g.,  harmonic  generations,  optical  frequency 
mixings,  and  optical  parametric  interactions)  usually  require  the  technique  of  phase  matching  in 
acentric  crystals  to  achieve  high  conversion  efficiency.  In  the  three-coupled-wave  interaction,  the 
conditions  of  energy  and  momentum  conservation  must  be  obeyed  are  written  as  follows; 

a>3  =  +  Wj  (1) 


where  0)  and  it  represent  angular  frequencies  and  wave  vectors,  respectively,  associated  with  the 
three  elecuomagnetic  waves.  Phase  matching  is  achieved  when  Ak  =0.  Eqn  (2)  can  be  expressed 
by  the  following  equation  using  refractive  indices  n  for  the  respective  frequencies: 

Qj/ijZj  =  +  Ojnji,  (3) 

— > 

where  •  ;•  the  unit  wmc:  parallel  ♦o  kj.  The  collinear  phase  matching  condition  requires  that  all 
three  wives  prcpaj.tte  Li  arallel,  i.e.,  tj  =  12  =  13.  which  reduces  eqn  (3)  to  the  form: 

0)3/13  =  0)2/12  ^ 

Eqn  (4),  written  in  terms  of  wavelength  becomes: 


Phase  matching  is  customarily  realized  by  using  the  birefringent  and  dispersive  properties  of 
anisotropiu  crystals.  Three  types  of  phase  matching  are  possible,  in  principle,  in  birefringent 
crystals  with  both  uniaxial  and  biaxial  characteristics.  They  are  Type  I,  where  O),  and  cOj  are  of 
the  same  polarization,  and  Types  II  and  III,  where  (O,  and  tOj  are  of  onhogonal  polarizations  (see 
Figure  1).  Types  II  and  III  phase  matching  differ  in  the  selection  of  co,.  In  Type  III  phase 
matching,  co,  is  smaller  than  coj.  As  such,  (Oj  will  have  larger  eigenpolarization  index  differences 
than  those  seen  by  co,.  In  general,  this  makes  Type  III  phase  matching  more  difficult  to  achieve. 
According  to  the  Biot-Fresnel  theorem  [5],  an  optical  wave  propagating  in  a  biaxial  crystal  with 
a  wave  vector  k  (0,4))  can  be  decomposed  into  two  components  with  mutually  perpendiculai 
polarization  directions  denoted  by  e,  (slow  ray)  and  Cj  (fast  ray),  respectively.  Beams  polarized 
in  these  and  only  these  directions  will  propagate  in  the  k  (0,4))  direction  without  change  in 
polarization  through  the  anisotropic  crystal.  But,  as  implied  above,  the  two  orthogonal 
polarizations  do  not,  in  general,  have  the  same  refractive  index,  n. 
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Figure  1.  Phase  Matching  Types  and  Associated  Polarizations 


In  uniaxial  and  biaxial 
crystals  there  exists  a 
direction  (one  in  the  case  of  a 
uniaxial  crystal  and  two  for  a 
biaxial  crystal)  in  which  the 
polarizations  do  have  equal 
refractive  indices.  This 
direction  is  called  the  optic 
axis  (see  Figure  2).  A 
uniaxial  crystal  is  essentially 
a  special  case  of  a  biaxial 
crystal  so  the  discusion  will 
continue  for  the  biaxial  case. 

In  biaxial  crystals,  the 
allowed  polarization 
directions  for  a  given 
direction  of  propagation  are 
as  follows.  The  spherical 
projection  in  Figure  2  is  used 
to  represent  the  directions 
inside  the  crystal.  Let  the 
direction  of  It  be  that  of  light 
propagation,  and  OA  and  OB 
the  two  optical  axes  in  the 
XZ  plane.  Spherical  triangle 
AITB  is  obtained  by  drawing 
the  great  circle  Alt  and  Bit 
on  the  sphere.  Then  the 
polarization  direction  of  the 
optical  wave  whose  velocity 
is  lower  than  the  other  is  that 
of  the  bisector  of  angle  AltB, 
i.e.,  the  direction  of  Cj.  The 
faster  ray  is  perpendicular  to 
e,,  i.e.,  the  direction  of  Cj. 

Hence  the  relation  n''  >  is 
satisfied  for  both  positive  and 
negative  biaxial  crystals, 
where  superscripts  e,  and  ej 
indicate  the  polarization 

directions  of  the  propagating  waves.  Note  that  these  polarization  directions  lie  on  the  surface  of 
the  sphere.  These  properties  are  usually  called  Biot-Fresnel’s  theorem.  (5,  6]. 

The  refractive-index  for  an  optical  wave  of  wavelength  A,  in  a  biaxial  crystal  satisfies: 
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Figure  2.  Spherical  Projection  Showing  the  Slow  (Cj)  and  Fast 
(62)  Eigenmodes  of  a  Wavefront  in  a  Biaxial  Crystal 
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(6) 


where  0  (the  angle  from  the  z-axis)  and  <j)  (the  angle  from  the  x-axis  in  the  xy  pla*'.'*’)  together 
define  the  internal  direction  of  the  incident  wave  normal;  riy,  and  n,  (the  principal  indices  of 
refraction)  are  found  using  Sellmeier  equations  (which  are  obtained  as  a  result  of  best  fits  to 
index  versus  X  measurements).  Note:  Figure  3  identifies  the  indicatrix  ellipsoid  for  biaxial 
crystals  in  which  the  crystal  indices  have  been  defined  as  n,  >  Rearranging  eqn  (6)  into 

a  form  more  easily  handled  and  making  the  following  substitutions  (to  increase  readability): 

-  *2  _  -2  L  -2  -2 

x-n^  ,  a  =  n^^,  6  =  ny„,  c 
=  sin^e  cos^cj),  =  sin^0  sin^4),  k^  =  cos^e 

we  obtain: 

(^  +  c)  -  ky  (a+c)  -  kf(a  +  Z>)jA;  +  (k^bc  +  kyOc  +  k^ab'j  =  0 

The  basic  form  of  eqn  (8)  is  x^  +  Bx  +  C  =  0.  Using  the  substitutions: 

B  =  -k^  (b  +  c)  -  ky  (a+c)  -  k^(fl  +  b) 

C  =  klbc  +  kyOC  +  klac 


(8) 


(9) 


and  reducing  results  in  the  following  form: 


In  eqn  (9),  the  terms  a,  b,  and  c  are  obtained 
from  Sellmeier  equations  with  the  appropriate 
Sellmeier  coefficients  for  a  given  wavelength 
(represented  by  “co”),  temperature  (we  used 
room  temperature)  and  crystal.  Eqn  (10)  is 
solved  for  each  of  the  three  wavelengths  used 
in  an  optical  parametric  oscillator  (OPO), 
optical  parametric  amplifier  (OPA)  or  second 
harmonic  generation  (SHG)  condition.  The  two 
real  solutions  are  the  index  of  refraction  for  the 
slow  and  fast  rays.  For  the  slow  ray  the  index 
of  refraction  is  represented  by: 


Figure  3.  Dielectric  Ellipsoid  for  a  Biaxial 
Crystal  with  Wavefront 
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(U) 


= 


n/2 


-  4C 


and  the  fast  ray’s  index  of  refraction  is  represented  by: 


"/  = 


(12) 


Of  the  six  available  wavelength  dependent  refractive  index  values  for  three-wave  interactions, 
only  three  are  suitable  for  determining  phase  match  conditions.  Table  1  indicates  the  appropriate 
slow  and  fast  ray  equation  combinations  which  satisfy  the  necessary  conservation  of  momentum 
condition.  Note  that  Table  1  has  a  or  a  after  the  phase  match  type.  This  indicates  a 
positive  or  negative  nonlinear  crystal  for  the  associated  phase  matching  type;  this  study  only 
deals  with  positive  nonlinear  crystals.  However,  using  Table  1  assignments  for  the  eigenmodes 
of  positive  and  negative  crystals,  the  MathCad®  phase  matching  worksheet  (see  Appendix  A) 
easily  handles  negative  nonlinear  crystals  without  modification. 


Table  1.  Potentially  Phase- Matchable  Combinations  of  Cj  and  Cj  Eigenmodes,  Where 


Polz 

SHG.SFM 

DFM32,  OPA32, 

DFM3,,  OPA31 

OPO 

Type 

In  Out 

In  Out 

In  Out 

In 

Out 

"KiXi  Xj 

Xj,X2  — >  Xy 

Xj.Xy  —4  X2 

Xj 

— >  Xy.Xj 

r 

61,61  62 

^2>^l 

62,61  ->  61 

^2 

— >  62,62 

r 

62,62 

62,62  ~^]€z 

^2*^2 

ir 

^2>^1  ^2 

62,62  6j 

^2 

-4  62,62 

n* 

62,62  62 

62,62  ^^62 

~4  62,62 

iir 

&l,^2  ^2 

^2,^2 

62,62  — ^  62 

«2 

-4  62,62 

in* 

62,62  -4  62 

—4  62,62 

NOTES: 


(1)  For  Birefringent  Positive  and  Negative  Crystals  in  the  Appropriate 
Nonlinear  Reporting  Frame  [4] 

(2)  Eigenmodes  e,  and  gj  (“e”  and  “o”  in  the  uniaxial  limit)  are  the  slow, 
fast  eigenmodes  for  positive  crystals  and  the  fast,  slow  eigenmodes  for 
negative  crystals.  Hence,  fast  and  slow  assignments  for  respective  waves 
i  =  1,  2,  3  are  identical  for  positive  and  negative  crystals,  and  X3  is  always 
the  fast  eigenmode. 
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Solving  the  appropriate  combinations  of  eqns  (11)  and  (12)  (identified  in  Table  1)  provides 
two  solutions  in  the  form  of  an  index  of  refraction  value  at  a  given  0,<})  angle  pair.  For  instance, 
for  Type  II  OPO  phase  matching  to  exist  (based  upon  the  IT  line  of  Table  1),  the  following 
equation  must  be  satisfied: 

Aj  ^2  Aj 

One  method  useful  in  identifying  the  index  of  refraction  equation  combinations  which  provide 
a  phase  matching  solution  is  the  graphical  method.  The  solutions  to  the  wavelength  dependent 
indices  (see  eqn  (10))  are  plotted  over  a  range  of  0  at  a  specific  (j).  The  initially  chosen  values 
of  4>  are  typicdly  the  stan  and  end  points  of  the  xy  plane  (i.e.,  <})  =  0°  or  90®).  The  plot  of  the 
resulting  solution  set  identifies  any  existing  phase  match  condition-represented  by  the 
meeting/crossing  point  of  ±e  two  plots  (i.e.,  if  a  phase  match  condition  exists  then  only  one  point 
in  the  solution  “set"  exists  for  which  a  combination  of  the  applied  wavelength’s  fast  and  slow 
eigenmoies  index  of  refraction  is  equal  to  the  index  of  refraction  for  the  generated  wavelength 
(defined  by  eqn  (13)).  This  method  is  explained  further  in  the  discussion  below. 


(13) 


I 


b.  Description  of  the  Computer-Based  Equations. 

See  Appendix  A  for  the  worksheet  listing.  ^^3 

MathCad*  4.0,  a  Microsoft*  Windows*  application 
from  MathSoft  Inc.,  was  used  to  develop  the 
computer-based  equations.  We  used  MathCad* 
because  it  required  a  minimum  of  programming, 
was  very  flexible,  and  the  worksheets  developed 
were  very  easy  to  debug.  MathCad*  has  proved  a 
very  powerful  tool  for  our  purposes. 

The  MathCad*  4.0  worksheet  restricts  the 
equations  to  the  first  quadrant  when  searching  for 
a  values  of  0.  The  resulting  solutions  are  then 
plotted  (i.e.,  the  phase  matched  combination  of  the 
/I; /Ay  and  n2lX2  terms  is  set  equal  to  the  associated 
/ij/Aj  curve).  The  resulting  phase  match  angle  pairs 
(0,(J))  may  then  be  used  in  the  calculation  of  the 
effective  nonlinear  coefficients  (d^^)  for  the  desired 
range(s)  and  the  associated  walkoff  angle(s).  Figure  4.  First  Octant  of  a  Proposed 

„  „  .  _  _  .  Positive  Nonlinear  Reporting  Frame. 

Sellmeier  Coefficients. 

The  worksheet  reads  the  Sellmeier  coefficients 
from  a  disk-based  ASCII  file.  This  file  is  composed  of  three  rows  of  four  columns.  Each  row 
represents  the  coefficients  for  the  x,  y,  and  z  axis,  respectively.  Each  colutxin  represents  the  A, 
B,  C,  and  D  Sellmeier  coefficients  for  the  following  form  of  the  Sellmeier  equation: 
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= 


B 


(14) 


A  + 


l-{Clkf 


-D-y} 

/ 


where  axis  represents  the  x,  y,  or  z  axis  (the  index  of  refraction  for  the  axes  or  «,)  at  the 
wavelength  specified  by  ‘A,’’.  MathCad  reads  the  file  of  Sellmeier  coefficients  (see  eqn  A-1)  into 
the  variable  S  (defined  as  S  :=  READPRN(^/ename,  where  filename  is  the  name  of  the  file 
containing  the  Sellmeier  coefficients  for  the  crystal  being  evaluated).  NOTE:  MathCad  reads  all 
data  files  from  the  same  subdirectory  to  which  the  worksheet  has  been  stored  (or  retrieved  from), 
unless  filename  is  an  “associated”  file  (see  MathCad  documentation). 

The  Identification  of  the  Three-Wavelengths. 

The  worksheet  expects  the  user  to  identify  two  of  the  three  wavelengths  from  which  MathCad 
calculates  the  third.  The  user  should  ensure  that  the  two  wavelengths  entered  result  in  MathCad 
calculating  a  positive  valued  wavelength.  See  eqns  A-3a  through  A-3c.  Next,  an  assignment  of 
the  three  wavelengths  as  X,,,  and  is  expected  (see  eqns  A-3d  through  A-3f)  so  that 
X,  >  X2  >  ^  results  of  the  assignments  are  displayed  immediately  below  each  assignment). 

Phase  Matching  Type  Category. 

Although  the  mathematics  are  essentially  the  same  when  calculating  phase  matching  for  SHG, 
OPA,  or  OPO  cases,  the  worksheet  (eqn  A-4)  provides  the  option  (for  those  who  so  desire)  to 
select  from  the  following  four  categories:  SHG/SFM,  OPA-32/DFM-32,  OPA-31/DFM-31,  or 
OPO.  The  category  selection  is  used  by  eqns  A-1  la  through  A-22f  to  identify  the  correct 
wavelength  dependent  index  of  refraction  equations  to  use  in  the  phase  match  solve  block  (see 
eqns  A- 13,  A- 14,  and  A-15)  for  Type  I,  H,  and  HI  phase  matching.  Note  however,  that  only  when 
it  is  time  to  calculate  the  nonlinear  coefficient,  d^g..  (see  appendix  B)  must  a  category  type 
actually  be  defined.  See  chapter  3  which  deals  with  the  calculation  of 

indices  of  Refraction. 

The  Sellmeier  coefficient  matrix,  the  variable  S  (see  eqn  A-1),  is  used  by  the  Sellmeier 
equations  (see  eqns  A-5a  through  A-5c)  to  calculate  tHe  indices  of  refraction  for  the  crystal  axes 
(/i,(X),  «,(X),  and  n,(X))  according  to  the  wavelengths  defined  in  eqns  A-3d  through  A-3f.  These 
values  are  stored  in  a  matrix  identified  by  the  variable  n  (see  eqn  A-5d).  This  variable,  n,  is 
saved  to  disk  in  eqn  A-5e;  WRITEPRN(/i/ename),  where  filename  is  the  name  of  the  file  to 
which  the  matrix  of  axes  indices  of  refraction  is  stored. 

Modifications.  If  the  Sellmeier  equations  provided  are  not  appropriate  (i.e.,  the  point  group  is 
different,  the  Sellmeier  equation  form  has  been  changed,  etc),  eqns  eqns  A-5a  through  A-5c  must 
be  modified  accordingly.  This  is  a  simple  edit  and  should  be  the  only  modification  required. 

Refractive-Index  Equation. 

The  solution  to  eqn  (6)  is  developed  with  eqns  A-6a  through  A-8f  The  solutions  of  the 
refractive-index  equations  are  given  by  eqns  A-9a  through  A-9f  and  represent  the  indices  for  the 
fast  and  slow  rays  for  each  of  the  three  wavelengths.  The  wavelength  dependent  versions  of  these 
equations  are  given  in  eqn  A-lOa  through  A-lOf.  They  are  then  arranged  by  eqn  A- 11a  through 
A-12o  into  formats  given  in  Table  1  for  the  various  SHG,  OPO,  and  OPA  combinations. 
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Phase  Match  Eauarions. 

The  first  step  in  finding  the  appropriate  phase  matching  equations  is  in  i.efining  the  wavelength 
dependent  refractive-index  equations.  These  are  arranged  so  that  two  ray  solutions  sets  are  set 
equal  to  a  third  solution  set  while  watching  the  fast  and  slow  ray  conditions  set  by  Table  1.  Eqns 
A- 11a  through  A-12o  define  the  wavelength  dependent  indices  while  equations  A- 13  through 
A- 15  define  the  appropriate  equations  to  use  in  the  phase  match  solve  blocks  as  defined  by  the 
category  (SHG,  OPA,  OPO,  etc.)  selected. 

The  phase  match  solutions  are  actually  found  by  the  solve  blocks.  Each  solve  block  (a  set  of 
equations  defining  a  solution)  is  composed  of  three  parts:  (1)  the  appropriate  phase  match 
equation  (see  eqns  A-  13a,  A- 14a,  and  A- 15a),  (2)  the  region  in  which  to  search  for  a  phase  match 
(this  worksheet  restricts  that  search  to  the  first  quadrant;  see  eqns  A- 13b,  A- 14b,  and  A- 15b),  and 
(3)  the  equation  solution  request  command  (i.e.,  FIND;  see  eqns  A- 13c,  A-14c,  and  A-15c). 

Phase  Match  Angle/Tvpe  Determination. 

Eqns  A- 16  through  A-22f  and  Figures  A-la  through  A-2c  are  used  obtain  and  display  the 
solutions  to  the  phase  matching  equations. 

To  begin,  we  first  determine  if  a  solution  set  exists  at  the  y-axis  in  the  xy  plane  (i.e„  <1>  =  90°). 
This  location  is  chosen  because,  since  n^>  tty  >  if  phase  matching  does  not  occur  in  the  yz 
plane,  it  will  not  occur  anywhere.  Eqn  A- 17a  sets  the  value  of  (})  to  90°,  eqn  A- 17b  sets  the  range 
of  0  (the  angular  range  from  the  z-axis  toward  the  y-axis)  over  which  to  search  for  a  phase  match 
angle,  and  eqn  A- 17c  convens  this  angular  range  to  radians  (MathCad  does  trigonometric 
calculations  in  radians).  Eqns  A- 18a  through  A-18f  assign  the  wavelength  dependent  refractive- 
index  equations  defined  in  eqns  A-1  la  through  A- Ilf  to  variables  for  use  in  Figures  A-la  through 
A-lc  to  graphically  display  the  results.  Type  I  phase  matching  equations  use  variables  A  and  B 
(see  eqns  A- 18a  and  A-18d,  respectively^  Type  n  phase  matching  uses  variables  C  and  D  (see 
eqns  A- 18b  and  A-18d,  respectively),  and  Type  HI  phase  matching  uses  variables  E  and  F  (see 
eqns  A- 18c  and  A-18f,  respectively).  When  one  line  in  any  of  the  figures  meets  or  crosses  the 
other  line,  a  phase  match  exists  at  that  point  for  the  associated  phase  match  type. 

The  value  of  the  variables  A  through  F  are  shown  above  each  figure.  These  displayed  values 
are  especially  useful  for  when  the  two  index  of  refraction  curves  appear  close  in  value,  but  an 
actual  phase  match  cannot  be  visually  determined.  There  should  be  at  least  six  digits  displayed 
to  the  right  of  the  decimal  point  to  ensure  correct  identification  of  a  phase  match  condition.  A 
phase  match  condition  exists  when  the  lower  number  (indicated  by  variable  B,  D,  and/or  F)  is 
equal  to  or  larger  than  the  upper  number  (indicated  by  variable  A,  C,  and/or  E). 

Eqn  A- 19a  through  A-22f  and  Figures  A-2a  through  A-2c  are  used  to  determine  where  phase 
matching  eacg  phase  match  type.  These  are  the  same  equations  and  figures  as  described  above, 
except  for  the  following: 

•  Eqns  A- 19a  through  A- 19c  identify  a  discrete  value  of  ({»  at  which  to  begin  a  phase 
match  search  (i.e.,  eqn  A- 19a  identifies  a  <j)  value  fr»r  Type  I,  eqn  A- 19b  identifies  a 
<{)  value  for  Type  n,  and  eqn  A- 19c  identifies  a  value  of  (j)  for  Type  III).  Eqn  A-21 
assigns  the  range  of  9,  defined  by  eqns  A-20a  through  A-2()c,  to  the  variable  0. 

•  The  values  for  eqns  A- 19a,  A-19b,  and  A- 19c  are  manually  varied  to  locate  the  start 
of  phase  match  region  for  each  phase  match  type  (for  those  phase  match  types  for 
which  a  phase  match  region  exists  at  (j)=90°  but  not  at  (})=0°). 
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The  start  of  the  phase  match  region  found  in  eqns  A- 19a  through  A- 19c  is  used  by 
eqns  A-23a,  A-24a,  and  A-25a  to  define  the  phase  matching  ranges  for  Types  I,  II,  and  ni, 
respectively.  If  a  particular  phase  match  type  does  not  exist,  the  value  of  0  for  all  values  of  (p  (0° 
to  90°)  are  automatically  zeroed  by  eqns  A-24b,  A-24c,  A-25b,  and  A-25c  (Types  II  and  in  only, 
since,  if  Type  I  does  not  exist  there  is,  by  definition  no  phase  matching  at  the  selected 
wavelengths).  NOTE:  MathCad  may  require  manually  zeroing  a  phase  match  region  (i.e.,  eqn 
A-24c  and/or  A-25c)  if  it  could  not  find  a  phase  match  for  that  phase  match  type. 

Saving  tlie  Phase  Match  Results. 

Once  phase  match  values  within  a  region  have  been  defined,  they  should  be  stored  to  disk  for 

(1)  easy  access  later  and/or  (2)  for  use  by  the  worksheet  which  calculates  the  and  walkoff 
angles.  This  is  performed  by  eqns  A-28a  through  A-29b.  To  aid  in  defining  the  data  file,  eqn 
A-27  identifies  the  phase  match  region  for  which  the  values  will  be  defined.  This  should  be  set 
to  a  range  of  0  to  90  to  allow  values  of  (p  over  the  entire  xy  plane.  Eqn  A-26  provides  conversion 
of  radians  to  degrees  used  by  eqns  A-28c,  A-28e,  A-28f,  A-28h,  and  A-28i  (the  worksheet  stores 
angle  data  to  file  in  degree  format,  not  radians).  Eqns  A-28a  through  A-28i  define  a  value  for 
each  range  index  value  (i.e.,  the  integer  value  of  (P).  These  values  are  then  arranged  into  a  data 
file  by  eqn  A-29a.  The  data  file  should  be  91  rows  (0-90,  as  defined  by  eqn  A-27)  of  7  columns 
(in  the  following  format:  X,,  Xj,  X3,  <p,  0(1),  0(11),  and  G(in),  where  the  first  3  values  are  the  3 
wavelengths  used,  die  next  value  is  the  integer  value  of  (p,  and  the  final  3  values  are  the  0  values 
for  Type  I,  II,  and  III  phase  matching  for  the  associated  value  of  (p).  The  variable  “data”,  eqn 
A-29a,  is  now  a  matrix  of  values  which  is  written  to  disk  by  eqn  A-29b.  This  equation  is  written 
as  WRl'l'EPRN  (filename)  :=  data,  'fjhs.xQ  filename  is  the  file  to  which  the  values  of  the  variable 
“data”  will  be  written.  Note  that  the  first  three  columns  of  data  each  have  identical  numbers.  This 
redundancy  is  required  to  allow  writing  the  data  as  a  single  matrix  with  information  which 
unambiguously  identifies  the  data  as  belonging  to  a  particular  wavelength  set. 

Final  Comments. 

(1)  The  value  of  <p  used  by  this  study  was  necessarily  restiicted  to  integer  values.  This  is 
fine  for  phase  match  regions  which  extend  across  the  entire  xy  plane  but  insufficient  for 
when  0°  <  <p  <  90°.  Therefore,  the  scratch  pad  region  provided  is  useful  for  determining 
the  value  of  <p  for  where  0  approaches  90°.  Specifically,  the  value  of  (p  is  varied  until  the 
effective  value  of  0  is  90°.  Note  that  the  equation  must  be  appropriate  to  the  phase  match 
type  (I,  n,  or  HI)  for  which  a  value  of  (p  is  desired  for  a  0=90°. 

(2)  MathCad  provides  a  method  of  converting  between  degrees  and  radians  automatically 
by  the  deg  function.  However,  since  the  same  function  is  used  to  conven  between  the  two, 
unless  carefully  used,  the  resulting  conversion  may  not  be  as  expected.  Therefore,  the 
worksheets  define  variables  which  perform  the  necessary  conversion  where  confusion 
would  most  likely  occur  when  using  the  MathCad  conversion  function. 

(3)  Phase  match  calculations  for  positive  or  negative  biaxial/uniaxial  crystals  are  identical 
with  the  supplied  worksheet.  No  modifications  are  necessary.  Not  until  c^culating  the  value 
of  d^ff  must  the  crystal  type  be  identified.  See  paragraph  3.b.ii  for  additional  information. 
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3.  EFFECTIVE  NONLINEAR  COEFnCIENT.  d,,f 


a.  Effective  Nonlinear  Coefficient  Theory. 

In  a  nonlinear  medium,  two  electric  fields  with  frequencies  co,  and  0)2  give  rise  to  induced 
polarizations  at  frequencies  O),  +  0)2  and  CD,  -  cOj.  The  energy  conversion  efficiency  of  a  three- 
wave  nonlinear  process,  such  as  SHG,  is  dependent  upon  the  size  of  this  induced  polarization, 
^e  second -order  polarization  P((D3),  induced  by  the  interaction  of  two  applied  electric  fields 
ECcOj)  and  ^(co,),  can  be  expressed  by  eqn  (15),  where:  4,-,  and  d*  are  the  unit  vector 


P(u,)  =  I  ie(®i)  I 

=  <i^|£(o,)|  |£(u,)| 

components  of  P((0,),  and  ^(O),),  respectively;  a,^  is  the  second-order  polarization  tensor 

of  the  nonlinear  crystal;  and  i,j,  and  k  each  represent  the  x,  y,  and  2  axes.  The  effective  nonlinear 
coefficient  is  a  function  of  the  angular  orientation,  the  crystal  symmetry,  and  the  values  of  the 
individual  components.  The  effective  nonlinear  coefficient  is,  thus,  given  by; 


In  a  three-wave  interaction  process,  is  equivalent  to  d,^,  so  that  it  can  be  expressed  by  a  3x6 
matrix  (as  d,).  For  a  more  complete  discussion  on  the  second  order  polarization  tensor  see 
paragraph  5.b. 

To  calculate  the  effective  nonlinear  coefficient  we  need  to  calculate  the  values  for  the  d„  d^, 
and  ^components.  The  value  of  dependent  upon  the  direction  of  the  H  and  since  we  know 
that  D  ±  K:  we  will  use  the  electric  displacement  vectors,  Ei^,(tOta)  to  calculate  d„  dy,  and  d*.  The 
electric  displacement  vectors,  ^,(0)a>).  where  r  =  1  (slow  rays)  or  r  =  2  (fast  rays)  for  these  two 
components,  and  their  projections  on  the  coordinate  axis  are  represented  by: 


where  i  =  x,  y,  and  z,  and  /n  =  Xj  and  Xj  [3).  The  slow  ray  (e,)  unit  vector  is  given  by: 

f  -COS0  cos4>  cos6^  +  sin<l) 


-COS0  sin<|>  cos6"  -  cos4>sin6” 


+sin0  cos6? 


/  Lin  \ 


(18) 
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and  the  fast  ray,  Cj,  unit  vector  given  by: 


m\ 


[^-cos0  cos4»  sin5*  -  sm<{>cosfi^ 
-COS0  sin<}>  sin5"  +  cos4>cos0P 
+sin0  sind” 


(19) 


where  57,  the  polarization  angle  (formed  between  the  two  eigenmodes,  e,  and  Cj,  and  the  zk 
plane)  is  expressed  as  a  function  of  0,  (j),  and  fit?  (see  eqn  (20)  as  well  as  Figure  2  and  Figure  4). 
See  references  [4]  and  [6]  for  further  details. 


5"  =  arctanj 


sin2<|>  cos0 


^cot^Q"  sm^0  +  sin^<|)  -  cos^0cos^<J) j 


*  0.5 


(20) 


The  angle  fit^,  given  by  eqn  (21)  below,  is  the  angle  formed  between  one  of  the  two  biaxial  optic 
axes  and  the  z-axis  (the  optic  axis  is  located  ±Q  from  the  z-axis)  in  the  xz  plane  [4].  The  unit 


Q*  =  tarcsin 


m 

(21) 


vectors  6?  (eqns  (18)  and  (19),  respectively)  denote  polarizations  of  electric  displacement  for 
biaxial  eigenmodes  e,  and  Cj  (see  Figure  2  and  Figure  4)  which  are  “e”  and  “o”,  respectively,  in 
the  uniaxial  limits.  In  a  principal-axis  system,  ^  and  0  are  related  as  shown  by  eqn  (22)  below. 


Vnl 

0 

0^ 

= 

0 

\fnl 

0 

0 

0 

l/< 

k 

(22) 
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Thus,  in  general,  the  unit  vectors  for  the  electric  fields’  two  polarization  components  are  given 

'»2(w«)J  ^23) 

Denoting  the  unit  vectors  of  ^/cdJ  as  d/co„  j  we  obtain  eqn  (24)  —  the  unit  vectors  for  the 
electric  field  polarization  component: 


d/oj  H  ar  = 


1 


Multiplying  these  unit  vectors  together  results  in  column  vectors  which  define  the  pseudo-vectors 
to  be  used  on  the  right  side  of  eqn  (16)  below).  Note:  m,  n,  and  o  represent  the  three  wavelengths 
(X„  X^,  and  X3)  involved  [3,  4]. 

Equation  (25)  defines  the  electric  field  polarization  column  vector  for  Type  I  SHG  phase 
matching  (rays  of  the  same  polarization  and  frequency): 


^  o’la."  = 


{<■4 

M 

K4 
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Equation  (26)  defines  the  electric  field  polarization  column  vector  for  Type  I  OPO  phase 
matching  (rays  of  the  same  polarization/different  frequency): 


„n  0 

o 


a.  y,'*  -nr® 

a"  y.®  a.  y,"  y.  ® 


Equation  (27)  defines  the  electric  field  polarization  column  vector  for  Type  Il/in  SHG  phase 
matching  (rays  of  different  polarization/same  frequency): 


^x,fl^x,t=2 
y*"  yr" 

“y.<=l“y.r=2 

y»**  yr" 


-  “f=l‘*r=2  m  m  m 

<t-l<f.2+<-l<».2, 

Equation  (28)  defines  the  electric  field  polarization  colunui  vector  for  Type  D/III  OPO  phase 
matching  (different  polarization  and  frequency): 


{^|.2(«o)J  =  <l<2  = 


y*"  y.® 

“j;f-l“a;»»2 

„n  y,® 
y."  y.® 

«  y.®  y®  y.® 

^.»=l‘*t.»=2  ■^^I.t=l^.r=2 
y_rt  _o  ,  _o 

®x,t-l^i,r»2  ^  ^z,f\^x,t=2 
„n  o  n  o 


From  these  results  we  are  able  to  obtain  the  general  expressions  for  the  effective  nonlinear 
coefficients  d,/I)  and  d,/II)  for  Types  I  and  II  phase  matching,  respectively,  using  the  definition 
for  d,ff  from  eqn  (16)  [3,  6].  Type  III  phase  matching  uses  the  d,/II)  equation.  Since  the 
difference  between  Type  11  and  III  phase  matching  is  that  the  polarizations  for  the  two  smaller 
frequencies  are  exchanged,  we  can  inspect  eqn  (27)  and  note  that  {d,^i(Gi„)}  fd,=2((0„)j  is  exactly 
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the  same  as  {d,^i((iii„)}  or  eqn  (28)  where  we  note  that  {d,^2(^o)}  is  exactly 

the  same  as  {d, ^2(^11)}  {d,:.i(oiJ}.  Therefore,  a  separate  equation  to  calculate  d,ff  for  Type  III  is 
unnecessary.  Note,  however,  though  the  form  is  the  same,  the  phase  matching  angles  for  Types 
II  and  ni  phase  matching  are  NOT  the  same  and  there  will  be  a  different  for  each  (e.g.,  see 
Figure  E-8  and  Figure  E-9). 


(29) 

(30) 

(31) 


=  a 


m  •  n  o 
Ut~2°iJk^J,t-l^k,f2 


(32) 


THe  subscripts  i,  j,  and  k  each  represent  the  x,  y,  and  z  axes.  It  should  be  noted  that  these 
equations  are  applicable  in  situations  where  ^  is  not  necessarily  perpendicular  to  1?. 


b.  Description  of  the  Computer>Based  Equations.  See  Appendix  B  for  the  worksheet 
listing. 


i.  Data  Retrieval/Initialization. 

This  worksheet  begins  by  reading  the  data  to  evaluate  from  three  ASCII  files.  These  files 
are  the  data  file  of  phase  match  information  (eqn  B-la),  the  2“*  order  polarization  tensor  (i.e.,  the 
djj  matrix;  eqn  B-lb),  and  the  indices  of  refraction  for  the  wavelengths  identified  by  the  first  three 
columns  of  the  data  file  (eqn  B-lc).  Note:  the  data  file  and  the  indices  of  refraction  file  must 
have  been  created  for  the  same  wavelength  values.  Next,  the  identification  of  the  origin  for 
matrices  is  identified  as  “0”  (i.e.,  row  1,  column  1  is  identified  as  0,0  while  row  3,  column  6, 
is  identified  as  2,5)  by  eqn  B-2. 


ii.  Positive/Negative  Nonlinear  Crystal  Identification. 

It  is  extremely  important  to  identify  whether  or  not  the  crystal  being  evaluated  is  positive 
or  negative.  Improper  identification  of  the  type  of  crystal  may  result  in  d^^  values  which  exactly 
opposite  that  expected  at  a  given  location  (i.e.,  a  maximum  rather  than  an  expected 
minimum).  Identifying  the  positive  or  negative  crystal  type  is  performed  by  eqn  B-3. 


iii.  Angle. 

Eqn  B-4  identifies  the  angle  <j)  for  which  a  phase  match  0  counterpart  is  desired.  The  angle 
pair  (0,<})),  when  found,  is  used  to  calculate  the  associated  and  walkoff  at  that  point. 
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iv.  Wavelengths. 

The  three  wavelengths  values  are  parsed  from  the  retrieved  data  file  and  assigned  to 
variable  A.,,  and  A.3  by  eqns  B-5a  through  B-5d.  Note  that  each  row  of  the  data  file  matrix 
begins  with  the  first  three  columns  being  the  wavelengths  at  which  the  data  was  collected. 

V.  Phase  Match  Angles. 

The  phase  match  angles  for  Type  I,  II  and  III  phase  matching  are  also  parsed  from  the  data 
file  by  eqns  B-6a  through  B-6c,  respectively.  For  proper  analysis  to  proceed,  the  worksheet  must 
know  which  phase  match  type  to  use  for  subsequent  calculations.  For  this  purpose,  the  value  of 
eqn  B-7a  is  manually  set  to  “1,”  “2,”  or  “3.”  Eqn  B-7b  then  takes  that  number  and  identifies  the 
associated  phase  match  type  matrix  (I,  II.  or  III).  This  identification  is  used  by  eqns  B-7c  and 
B-7e  to  assign  values  for  0  and  (|),  respectively.  These  assigned  values  are  redefined  specifically 
as  0  and  <J)  by  eqns  B-7d  and  B-7f,  respectively.  These  are  the  assignments  used  by  the  equations 
which  follow.  Eqn  B-8  selects  the  desired  category  type  (SHG/SFM,  OPA/DFM  or  OPO)  for 
which  calculations  are  desired.  Note  that  if  the  first  two  wavelengths  are  equal  the  worksheet 
automatically  selects  category  1  (i.e.,  SHG)  regardless  of  the  category  manually  selected. 

vi.  Miscellaneous. 

Eqns  B-9  through  B-16  are  discussed  in  the  theory  (see  chapters  2  and  3)  above.  In  short, 
however,  the  following  is  calculated: 

(1)  Eqn  B-9  calculates  the  value  of  Q  (the  angle  from  the  z-axis  to  the  optic  axis). 

(2)  Eqn  B-lOa  is  used  by  eqn  B-lOc  which  calculates  the  value  of  5  (the  polarization 
angle  formed  by  the  direction  of  the  slow  ray,  e,,  with  respect  to  the  surface  of  the  plane  formed 
by  the  wavefront  vector,  k  ,  and  the  z-axis).  We’ve  arbitrarily  restricted  the  value  of  15,  |  to 
^  90°.  The  sign  of  5  is  the  same  as  that  given  by  the  numerator  in  eqn  B-lOa.  Therefore,  when 
the  value  of  the  numerator  in  eqn- 10a  is  less  than  zero,  a  value  of  n  is  added  to  the  value  of 
atan(tan(28)  before  halving  it  to  get  the  correct  5,.  This  is  done  by  eqn  B-lOb. 

(3)  Eqns  B-lla  through  B-lll  are  setup  for  Eqns  B-llm  and  B-lln  which  calculate 
the  electric  displacement  unit  vectors  for  the  slow  rays  and  the  fast  rays,  respectively. 

(4)  Eqn  B-12  calculates  the  electric  field  polarization  components. 

(5)  Eqns  B-13a  and  B-13b  calculate  the  unit  vectors  for  the  polarization  components 
of  the  electric  fields  for  the  slow  r.  .ys  and  fast  rays,  respectively. 

(6)  Eqns  B-14a  through  B-14c  identify  the  proper  columns  of  eqns  B-13a  and  B-13b 
to  use  when  calculating  the  values  of  a,  (polarization  unit  vector  for  the  “generated”  wavelength), 
a^,  and  a^  (polarization  unit  vectors  for  the  “applied”  wavelengths). 

(7)  Eqns  B-15a  and  B-15d  identify  the  polarization  unit  vectors  a;,  a^,  and  a,j. 

(8)  Eqns  B-16a  and  B-16b  identify  the  fast  and  slow  column  vectors,  a^a^  (a.k.a., 
pseudo- vectors).  Eqn  B-16c  selects  the  appropriate  column  vector,  based  upon  the  phase  match 
type  selected,  to  multiply  with  the  second  order  polarization  tensor,  d,j. 

vii.  Calculating  the  Effective  Nonlinear  Coefficient,  d.^. 

( 1 )  The  calculation  for  SHG/SFM  and  OPO  d^g  results  in  exactly  the  value.  Also,  the 
calculation  for  d.^  for  OPA-32  results  in  the  same  value  as  for  OPA-31.  However,  if  the  same 
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wavelengths  used  for  SHG/SFM  or  OPO  are  used  to  calculate  for  an  OPA  case,  the  will 
be  different.  Therefore,  it  is  imponant  to  ensure  that  the  proper  category  is  selected. 

(2)  Eqn  B-17  calculates  phase  match  type  (i.e.,  I,  II,  or  III)  d^jj. 
viii.  Saving  the  Resulting  Calculations.  The  data  associated  with  the  d^ff  and  walkoff 
angle  calculations  are  compiled  into  a  single  row  matrix.  This  matrix  is  composed  of  the  three 
wavelengths  used,  the  phase  match  angles  (0,<|)),  the  value,  and  the  walkoff  angles  (for  the 
discussion  on  walkoff  angles  see  chapter  4)  for  the  phase  match  type  for  which  this  data  was 
calculated. 

Eqn  B-19a  provides  a  variable  for  converting  from  radians  to  degrees.  Eqns  B-19b  through 
B-I9d  selects  the  correct  fast  or  slow  walkoff  angle  for  wavelengths  A,,,  and  A.3.  The  matrix 
is  assigned  to  the  variable  “data”  (eqn  B-20a)  which  can  then  be  used  to  append  to  an  existing 
matrix  file  (as  an  additional  row  of  data).  Tliis  equation  (eqn  B-20b),  APPENDPRN(^/cname) 
:=  data,  defines  filename  as  the  file  to  which  to  append  the  data,  and  data  is  the  variable 
containing  the  single  row  matrix  of  data  described  above. 


16 


WALKOFF  ANGLE 


a.  Theory  for  the  Walkoff  Angle  Calculations. 

The  walkoff  angle  of  a  light  wave  in  a  nonlinear  crystal  is  the  angle  between  the  wave  vector 
It  and  the  Poynting  vector  » .  It  is,  therefore,  also  the  angle  between  the  electric  field  ^  and  the 
electric  displacement  vector  0  [3].  See  Figure  5  below.  To  illustrate,  consider  the  following 
which  says  that  the  magnitude  of  the  electric  field  is  the  same  as  the  magnitude  of  the  electric 
displacement  vector  times  the  electric  field  polarization. 


where: 


\E\  =  P\D\ 


where  i  =  jc,  y,  and  z;  b,  is  a  unit  vector  component  of  the 
electric  displacement  vector  for  the  slow  and  fast  rays  (see 
eqns  (18)  and  (19));  and  /i,-  is  index  of  refraction  for  each 
crystal  axis  for  a  given  material.  Since  the  following 
relationship  is  true  by  definition  of  the  dot  product. 


cosp  = 


M'S 

\E\\D\ 


Figure  S.  Relationship  Between 
ic",  f,  and  ^ 


and  the  fact  that  ^  1 and  S  iJi,  finding  the  angle  between  the  wave  vector  and  the  Poynting 
vector  is  the  same  as  finding  the  angle  between  the  electric  field  and  the  electric  displacement 
vector.  Now,  let  Ej  =  |El  a;.  And,  using  eqns  (17)  and  (23),  we  can  say: 

=  j:(|£|a,X|DlJ.,)  (36) 


which  leads  to: 


E-3  =  |£|  |D|X; 

I  ''  n, 

D  ^  2 


and  therefore  we  set: 
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Now,  using  eqn  (35),  we  reduce  eqn  (38)  to  a  trigonometric  equation: 


cos  p 


(38) 


(39) 


which  reduces  to: 


which  is  the  form  used  in  eqn  (41)  (the  form  used  in  the  MathCad  worksheet)  below. 
Consequently,  the  walkoff  angles  of  the  light  waves  polarizing  in  e,,  (where  /  =  1  or  r  =  2)  in  a 
nonlinear  crystal  can  generally  be  expressed  as  [3]: 


b.  Use  of  the  Computer>Based  Equations. 

See  Appendix  B  for  the  listing.  Eqn  B-23  calculates  the  walkoff  angles  for  the  slow  and  fast 
rays.  The  calculation  is  straight  forward  and  requires  no  explanation  beyond  that  described  in  the 
theory  above.  There  are  no  user  modifications,  additions,  or  changes  required.  However,  as  in  the 
case  of  the  calculation  of  the  proper  data  files  must  be  retrieved  and  the  phase  match  type 
(I,  n,  or  III)  must  be  identified.  See  paragraph  3.b  for  more  information  on  data  retrieval  and 
initialization  within  MathCad. 
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5.  LINEAR  AND  NONLINEAR  PROPERTIES  OF  KTP  ISOMORPHS. 


a.  List  of  Properties  Used  in  Calculations. 

Table  2  and  Table  3  summarize  the  Sellmeier  coefficients  and  nonlinear  coefficients  of  the 
KTP  isomorphs,  respectively.  All  calculations  in  the  phase  matching  and  d^^^/walLoff  MathCad 
worksheets  use  the  information  supplied  by  these  tables.  Phase  matching  curves  derived  fc»  these 
isomorphs  are  qualitative  for  wavelengths  longer  than  1.15  pm  because  of  less  experimental 
testing. 


Table  2.  Sellmeier  Coefficients  of  KTP  Isomorphs 


Axis 

KTP 

RTP 

KTA 

RTA 

CTA 

A, 

2.11460 

2.15559 

2.11055 

2.22681 

2.34498 

B, 

0.89188 

0.93307 

1.03177 

0.99616 

1.04863 

c. 

0.20861 

0.20994 

0.21088 

0.21423 

0.22044 

D, 

0.01320 

0.1452 

0.01064 

0.01369 

0.01483 

A, 

2.15180 

2.38494 

2.38888 

1.97756 

2.74440 

B, 

0.87862 

0.73603 

0.77900 

1.25726 

0.70733 

c, 

0.21801 

0.23891 

0.23784 

0.20448 

0.26033 

P. 

0.01327 

0.01583 

0.01501 

0.00865 

0.01526 

A. 

2.31360 

2.27723 

2.34723 

2.28779 

2.53666 

B. 

1.00012 

1.11030 

1.10111 

1.20629 

1.10600 

c. 

0.23831 

0.23454 

0.24016 

0.23484 

0.24988 

Dx 

0.01679 

0.01995 

0.01739 

0.01583 

0.01711 

Sellmeier  coefficients  of  best  fit  to  index  data,  with 

n]  =  A,  +  (B,  I[1  -  (C,  iXfj)  -  D,‘,  where  i  =  x,  y,  and  z,  respectively 

(corresponding  to  the  crystallographic  axes  a,  b  and  c.)  [7} 
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Table  3.  Nonlinear  Optical  Properties  of  KTP  Isomorphs 


d, 

KTP 

RTP^ 

KTA 

RTA^ 

CTA^ 

dn 

16.9 

17.1 

16.2 

15.8 

18.1 

du 

4.4 

4.1 

4.2 

3.8 

3.4 

dji 

2.5 

3.3 

2.8 

2.3 

2.1 

d24 

3.6 

3.6 

1.7  X  dii  of 
KTP(a  1.3pm 

3.6 

3.6 

d,s 

1.9 

1.9 

1.2  X  d,s  of 
KTP@  1.3pm 

1.9 

1.9 

SHG  d,j  coefficients  at  1.064  fim  (in  pmlV;  experimental  uncertainties:  ±  10%  for  djj, 
d24,  and  dn,,  ±  20%  for  d32  cutd  d^  [7} 

f  d,s  and  d24  were  not  measured  by  Cheng,  Cheng,  and  Bierlein  (see  reference  [7]) 

—  so  the  KTP  value  was  used  here.  However,  they  did  state  that  within  experimental 
errors,  there  is  no  detectable  variation  of  the  d  coefficients  among  the  isomorphs 

b.  Nonlinear  (dy^)  Coefficients, 
i.  Basic  Tensor  Background. 

Eqn  (15),  shown  again  as  eqn  (42),  represents  the  nonlinear  (2nd  order)  polarization.  Since 
no  physical  significance  can  be  attached  to  an  exchange  of  E„,  and  E„2  it  follows  that 

PCo,,)  =<J^|£(Ul)||B(Uj)|  («) 

This  means  we  can  replace  the  subscripts  kj  by  a  single  symbol  (e.g.,  kj  may  be  represented  by 
only  j  term  without  any  loss  of  information  or  accuracy)  according  to  the  contraction: 

xc  =  1  yy  =  2  2z  =  3 
yz  =  zy  -  4  jy-  =  zx  =  5  .xy  =  yx  =  6 

The  resulting  d,^  tensor  forms  a  3  x  6  matrix  (down  from  a  multidimensional  3x3x3  array  that 
operates  on  the  E*  column  tensor  to  yield  P  according  to  eqn  (43).  The  contracted  d,^  tensor  obeys 
the  same  symmetiy  restrictions  as  the  d,^^  tensor,  and  in  crystals  of  a  given  point-group  symmetry 
it  has  the  same  form.  [9] 
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£* 

^12  ^15  ^ 

^  *43  *44  ^  ^  * 

<4l  ^2  ^3  ^  4jS  <46 

2£,£, 

2£,£, 

ii.  Tensor  Matrix  for  KTP  Isomorphs. 

The  KTP  family  of  crystals  belongs  to  a 
noncentrosymmetrical  class  called  orthortiombic.  In  that 
class,  the  nonlinear  optical  tensor  is  of  class  ntm2.  It  has 
the  form  shown  in  Figure  6.  The  symbol  indicates  a 
zero  modulus  and  indicates  a  nonzero  modulus.  Notice 
that  according  to  the  given  form,  KTP  has  five  nonzero  d 
coefficients;  djj,  d22,  d^,,  and  djy  The  coefficient  djj  is 
approximately  9  times  greater  in  magnitude  than  the 
diagonal  coefficients  {d,^  and  whose  values  are  within  Figure  6.  An  mml  Crystal  Class 
a  factor  of  two  of  each  other.  The  effective  nonlinearity  Tensor  Arrangement, 
for  Type  II  (or  Type  HI)  in  the  principal  ttianes  is  given  by 
Table  4  below  (note:  for  this  to  be  true,  E  must  be  parallel  to  1^). 

Table  4.  Effective  Nonlinearity  for  Type  II  SHG  in  the  Principal  Planes 
PLANE  d^ff 

x-z  d24  sin9 

y-z  d,5  sinG 

x-y  d/5  sin^(|)  +  d24  cos^tj) 

Thus,  the  larger  djj  coefficient  has  only  a  very  minor  influence  on  the  magnitude  of  d,^.  Also, 
since  ^24  is  greater  than  d/5,  Type  IIAII  phase  matching  is  generally  more  efficient  in  the  jcz  plane 
than  in  the  yz  plane  [7]. 

c.  Caveats. 

No  study  is  complete  without  identifying  baselines  used  to  obtain  or  calculate  data.  A  number 
of  assumptions  were  made  in  this  study  because  of  the  newness  of  the  crystals  being  studied  and 
because  of  the  controversy  surrounding  certain  values  of  those  crystals. 
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i.  (i-Coefficients. 

The  isormorph  rf-values  used  are  the  most  recent  available  at  this  writing  (see  Reference 
f7]).  Only  the  values  for  KTP  and  KTA  are  well-known,  although  still  apparently  not  fully  agreed 
upon.  These  values,  supported  by  [10],  result  in  a  which  is  roughly  half  of  earlier  reported 
values.  The  di5  and  d24  values  for  RTP,  RTA,  and  CTA  were  not  known  (or  published)  at  the 
time  of  this  writing  so  the  associated  values  for  KTP  were  substituted.  For  this  reason,  the  results 
calculated  for  those  crystals  may  slightly  differ  from  their  true  values.  However,  according  to  [7], 
within  experimental  error  there  is  no  detectable  variation  of  the  d  coefficients  among  the 
isomorphs.  Further  study  is  required. 

ii.  Miller’s  Rule. 

Miller  showed  that  although  nonlinear  coefficients  may  vary  considerably  among  nonlinear 
crystals,  a  parameter  (A,^  relating  the  nonlinear  and  linear  optical  susceptibilities  is  constant  to 
within  an  order  of  magnitude.  This  constant  is  defined  by  eqn  (44)  [11,  12]: 


€o  ["/(Wj)  -  1]  [«/(w2>  ■  -  1] 


(44) 


The  values  of  /j„  /jj,  and  represent  the  indices  of  refraction  along  the  principal  axes. 
Experimental  values  for  d,y*  used  in  this  study  were  at  the  “base”  3-wave  interaction  of  1.064  pm 
SHG  [7].  Adjustments  to  the  d-values  for  use  at  other  wavelengths  are  made  by  application  of 
Miller’s  Rule.  The  following  background  information  should  help  considerably  in  properly 
applying  Miller’s  Rule. 

As  stated  in  paragraph  5.b  above,  the  d,j  tensor  is  a  contracted  form  of  a  multidimensional  set 
of  3x3x3  matrices  (dy*).  The  resulting  form  (identified  in  Figure  6)  is  displayed  in  eqn  (45)  below 
with  the  tensor  matrix  identification.  This  form  identifies  the  correct  indices  of  refraction  to  use 


xxx  xyy  xzz  xyz  x*z  xxy 
yxx  yyy  yzz  m  yxz  yxy 
fxx  xyy  zzz  zyz.  zxz  zxy, 


(45) 


when  applying  Miller’s  Rule.  Therefore,  the  five  nonzero  d-values  for  KTP  and  its  isomorphs  are 
more  correctly  identified  as:  d/.,  =  d^,  d^  =  d^^.  dj,  =  d^,  dj,  =  d^,  and  d^j  =  d.„.  In  applying 
this  information,  note  that  the  first  value  indicates  that  the  refractive  index  value  of  the  identified 
axis  at  (Oj  (the  calculated  wavelength,  i.e.,  (Oj  =  to,  +  to,),  the  second  value  is  for  the  refractive 
index  value  of  the  identified  axis  at  (O,,  and  the  same  thing  applies  for  co,.  Also,  note  that  the 
indices  of  refraction  at  a  specific  wavelength  are  applied  according  to  0)_,  <  co_i  <  to,.  Each  letter 
in  the  tensor  matrix  values  identifies  the  axis  of  the  appropriate  index  of  refraction.  NOTE: 
because  d.^^  =  d,*^,  the  last  two  indices  of  refraction  may  be  exchanged  without  changing  the  value 
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of  dij  (c.g.,  dis  =  d^  or  d,s-d^).  For  example,  in  applying  Miller’s  Rule  to  calculate  for  d,s, 
eqn  (44)  would  appear  as  follows: 


= 


‘15 


eo*(«*(“3)  -  -  l) 


(46) 


To  obtain  the  corrected  ^/-value  at  the  new  wavelengths,  eqn  (44)  is  rearranged  as: 


"  1H«/(«2)  -  1][»4(Mi)  -  1] 


(47) 


Taking  Miller’s  Rule  into  account  for  SHG  in  KTP  of  1.064  pm  fundamental  input  (the  three- 
wave  interaction  for  which  the  d-values  of  KTP  are  given)  and  making  adjustments  to  d-values 
at  the  SHG  and  OPO  wavelengths  used  in  this  study,  we  obtained  tne  data  presented  in  Table  5 
below. 

Miller’s  Rule  corrections  to  (1.064/1.064/0.532  pm)  d-values  for  the  other  interactions 
considered  represent  changes  of  22%  or  less.  However,  since  a  considerable  amount  controversy 
exists  as  to  exactly  what  are  the  true  values  of  d  and  experimental  measurements  are  only  good 
to  10-20%,  the  importance  of  Miller’s  Rule  is  considerably  minimized.  Therefore,  Miller’s  Rule 
corrections  were  not  used  in  this  study. 


Table  5.  Miller’s  Rule  Applied  to  1.064  pm  KTP  d  Values 


KTP 

Wavelengths  (pm) 

(/•value 

Uncertainty 

Range 

(pm/V) 

Experimental 

Uncertainty 

171  1 

1.064/1.064/ 

0.532 

13SH3SI 

0.675 

4.043/1.444/ 

1.064 

3.17/1.60/ 

1.064 

d  (pn 

i/V) 

16.9 

15.72 

13.25 

13.87 

15.2-18.7 

WESMi 

Wm 

4.4 

4.11 

3.51 

3.66 

4.0-4.8 

BM 

2.5 

2.34 

2.00 

2.08 

2.3-2.8 

hh 

3.6 

3.39 

2.88 

3.01 

2.9-4.3 

±20% 

19 

1.9 

1.79 

1.53 

1.60 

1. 5-2.3 

±20% 

iii.  Sellmeier  Coefficients. 

There  is  a  difference  in  the  Sellmeier  coefficients  depending  upon  whether  KTP  was  flux 
or  hydrothermally  grown.  This  fact  could  have  a  small  impact  upon  the  phase  matching  results. 
The  same  situation  might  also  apply  to  the  KTP  isomorphs. 
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6.  ISOMORPH  COMPARISONS 


a.  Introduction. 

The  interacting  frequencies  (wavelengths)  considered  for  this  study  were  of  interest  for 
on-going  in-house  experiments.  Second  harmonic  (SHG)  calculations  used  1.064  and  1.35  pm  as 
the  fundamental  wavelengths  since  these  are  commonly  available  Nd-laser  wavelengUis.  Optical 
parametric  oscillator  (OPO)  calculations  considered  1.064  pm  as  the  pump  wavelength,  1.444  and 
1.60  pm  as  the  signal  wavelengths,  and  4.043  and  3.17  pm  as  the  idler  wavelengths,  respectively. 
These  wavelengths  were  chosen  because  of  current  investigations  using  OPO  and  DFM  to 
generate  mid-infrared  wavelength  laser  sources. 

b.  Phase  Matching  Characteristics. 

The  range  over  which  phase  matching  is  possible  for  the  four  three-wave  interactions  is  shown 
in  tabular  form  in  Appendix  C.  Plots  of  phase  matching  are  shown  in  Figure  7  through  Figure  15. 
Calculations  were  made  at  1°  intervals  of  ({>;  the  symbols  are  shown  every  10°  for  labeling.  All 
plots  show  a  decrease  in  0  with  an  increasing  (j).  The  overall  change  in  0  over  the  available  phase 
match  range  is  given  by  the  value  F.  Clearly  a  value  of  F  =  0  indicates  a  uniaxial  case.  Thus, 
this  value  is  an  indicator  of  how  close  to  uniaxial  each  biaxial  crystal  is  over  a  given  phase 
match  region.  Table  6  through  Table  10  summarize  the  angular  phase  matching  ranges. 

i.  Second  Harmonic  Generation. 

(1)  Type  I.  Figure  7  and  Figure  8  show  Type  I  SHG  phase  matching  for  the  KTP 
isomorphs.  Type  I  phase  matching  is  possible  with  all  isomorphs  throughout  the  ^-angle  range. 
RTP  and  CTA  have  much  greater  F’s  than  the  other  isomorphs  by  approximately  factors  of  3  and 
2.5,  respectively. 


Table  6.  Type  I  SHG  Variations  in  0 


Cfystal 

1.064/1.064^.532 

1.35/1.35/0.675 

(Typ«l 

SHG) 

Rangt  of  6 

r 

Rangt  of  6 

r 

KTP 

49.41*  •  39.20° 

0.11 

40.25°  -  29.57° 

0.12 

RTP 

55.83°  -  32.17° 

0.30 

47.12°- 18.82° 

0.31 

KTA 

51.47°  -  44.69° 

0.08 

40.76°  -  33.79° 

0.08 

RTA 

57.66  -  47.48° 

0.11 

45.08  -  35.06° 

0.11 

CTA 

75.50°  -  52.22° 

0.25 

56.37°  -  34.26° 

0.25 
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<t>  (Deg) 

Figure  7.  SHG  Phase  Matching,  Type  I, 
1.064/1.064/0.532  pm. 


0  10  20  30  40  50  60  70  80  00 
<p  (Deg) 

Figure  8.  SHG  Phase  Matching,  Type  I, 
1.35/1.35/0.675  pm. 


(2)  Type  II.  Figure  9  and  Figure  10  show  Type  11  SHG  phase  matching  curves.  None 
of  the  three  arsenate  isomorphs  exhibited  Type  II  phase  matching  for  SHG  using  1.064  pm  at  any 
orientation.  Even  the  phosphates  had  a  limited  range  of  ^  angles  at  1.064  pm.  At  1.35  pm,  SHG 
Type  II  phase  matching  was  more  easily  achieved  with  all  five  isomorphs  phase  matchable  and 
only  CTA  having  a  restricted  <()  range. 


Table  7.  Type  fl  SHG  Variation  in  0 


jjRQH 

1.064/1.064/0.532 

1.35/1.35/0.675 

mSM 

Rang*  of  e 

r 

Rang*  of  e 

r 

KTP 

gO-OO”  -  68.67'’ 

0.33 

49.11”  -  39.42” 

0.11 

RTP 

90.00”  -  64.62” 

0.53 

65.95”  -  43.91” 

0.24 

KTA 

67.67”  -  54.67” 

0.08 

RTA 

— 

70.40  -  59.73” 

0.12 

CTA 

90.00”  -  72.04” 

0.49 

(3)  Type  III.  For  SHG  phase  matching  purposes,  this  type  is  exactly  the  same  as  Type 
II  since  the  two  input  photons  have  the  same  wavelength.  As  soon  as  the  degeneracy  is  lifted 
though,  they  become  different  and  this  will  become  significant  for  the  OPO  calculations  discussed 
below. 
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f  ''iY  '■* '•■"  '  '  .  '  ^  ‘  'X*** 


0  (Deg) 

Figure  9.  SHG  Phase  Matching,  Type  Il/ni, 
1.064/1.064/0.532  pm. 


0  (Deg) 

Figure  10.  SHG  Phase  Matching,  Type  Il/ni, 
1.35/1.35/0.675  pm. 


ii.  Optical  Parametric  Oscillator. 

(1)  Type  I.  Figure  1 1  shows  Type  I  phase  matching  for  an  OPO  pumped  at  1.064  pm 
with  1.444  pm  signal  and  4.043  pm  idler  wavelengths.  Figure  12  shows  the  phase  matching 
results  for  the  1.60  pm  signal  and  3.17  pm  case.  Like  the  SHG  results.  Type  I  phase  matching 
is  possible  for  all  isomorphs  at  all  ^  angles. 


Table  8.  Type  I  OPO  Variation  in  6 


Crystal 

4.043/1.444^.064  pm 

3.17/1.60/1.064  pm 

(Typ«l 

OPO) 

Rang#  of  6 

r 

Rang*  of  e 

r 

KTP 

39.24”  -  29.06” 

0.11 

37.62” -27.03” 

0.12 

RTP 

47.61”  - 19.74” 

0.32 

45.85”  - 15.77” 

0.33 

KTA 

37.30”  -  32.98” 

0.05 

35.95”  -  30.92” 

0.06 

RTA 

41.94”  -  27.80“ 

0.16 

40.35”  -  26.65” 

0.15 

CTA 

53.80”  -  29.52” 

0.28 

51.89”  -  26.41” 

0.24 
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0  (Deg) 

Figure  11.  OPO  Phase  Matching,  Type  I, 
4.043/1.444/1.064  pm. 


^  (Deg) 

Figure  12.  OPO  Phase  Matching,  Type  I, 
3.17/1.60/1.064  pm. 


(2)  Type  H.  Type  II  OPO  phase  matching  curves  are  shown  in  Figure  13  and 
Figure  14.  Like  Type  I,  Type  II  phase  matching  is  achievable  over  all  (])  values  in  all  isomorphs. 
The  basic  difference  is  that  the  Type  II  0  phase  matching  angles  vs.  (|)  are  a  few  degrees  higher 
than  the  corresponding  0’s  for  Type  I. 


Table  9.  Type  II  OPO  Variation  in  0 


Crystal 

4.043/1.444^.064  pm 

3.17/1.60/1.064  pm 

(TypsM 

OPO) 

Rang*  of  6 

r 

Rang*  Of  8 

Hi 

KTP 

44.33"  -  35.47" 

0.10 

44.85"  •  35.96" 

0.10 

RTP 

52.79"  -  28.08" 

0.27 

53.16"  •  28.67" 

0.27 

KTA 

45.72"  -  34.82" 

0.12 

46.01"  •  36.30" 

0.11 

RTA 

42.48"  -  39.65" 

0.03 

44.43" -40.17" 

0.05 

CTA 

61.59"  -  40.06" 

024 

6^60"  -  41.18" 

024 

<f>  (Deg) 

Figure  13.  OPO  Phase  Matching,  Type  II, 
4.043/1.444/1.064  pm. 


<j>  (Deg) 

Figure  14.  OPO  Phase  Matching,  Type  II, 
3.17/1.60/1.064  pm. 
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(3)  Type  III.  No  Type  HI  phase  matching  could  be  achieved  by  any  isomorph  at 
4.034/1.444/1.064  pm.  Because  dispersion  decreases  with  increasing  wavelength. 


(48) 


Adding: 


(49) 


yields: 

^lUrwi^bOtr)  ^  ^jaM^UOtt)  *  ^slow{^tltitat) 


(50) 


Here,  the  left  side  is  the  index  sum  used  for  Type  II  and  the  right  side  is  the  index  sum  for  Type 
III.  Usually,  lack  of  phase  matching  means  that  the  sum  is  insufficient  to  match  the  pump  index. 
Thus,  Type  II  phase  matching  is  “easier”  since  it  always  has  a  larger  sum  and,  for  the  above 
case,  (4.034/1.444/1.064  pm)  Type  II  works  while  Type  III  does  not. 

For  the  3.17/1.60/1.064  pm  case.  Type  ID  phase  matching  was  possible  for  all  isomorphs 
except  for  CTA.  Since  CTA  has  the  largest  6  value  at  a  given  (])  for  all  types  of  phase  matching, 
it  logically  would  be  the  first  to  fail  to  phase  match  for  Type  III  where  higher  0  values  are 
needed.  Only  KTP  could  be  phase  matched  for  all  <})  values. 


Table  10.  Type  QI  OPO  Variation  in  0 


Crystal 

(Typslll 

OPO) 

3.17/1.60/1.064  pm 

Bangs  of  e 

r 

KTP 

77.36'  -  65.90° 

0.13 

RTP 

90.00'  -  66.88' 

0.52 

KTA 

90.00' -  71.37' 

0.25 

RTA 

90.00' -  81.52' 

0.27 

CTA 

iii.  Comments. 

The  largest  variations  in  0  were  present 
for  the  crystals  RTF  and  CTA  while  KTA, 

KTP,  and  RTA  showed  much  smaller  0 
variations.  Thus,  KTA,  KTP,  and  RTA  are  very 

uniaxial-like  in  nature  (i.e..  =  n^)  and  the  RTP  and  CTA  crystals  are  more  biaxial  in  nature 

(i.e.,  /I, «  nj. 


(p  (Deg) 

Figure  15.  OPO  Phase  Matching,  Type  III, 
3.17/1.60/1.064  pm. 
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c.  Effective  Nonlinear  Coefficients 

Appendix  E  shows  calculations  for  all  SHG  and  OPO  interaction  types  considered.  Each 
plot  shows  dfff  as  a  magnitude  above  the  0(t)  phase  matching  plane.  All  plots  in  Appendix  E  use 
the  same  vertical  scale  for  ease  of  comparison.  The  maximum  for  each  isomorph  at  the  SHG 
and  OPO  wavelengths  according  to  phase  matching  type  is  shown  in  Table  1 1.  Recommendations 
for  any  particular  crystal  in  this  section  are  related  only  to  the  largest  value  of  d^g  found. 

i.  Second  Harmonic  Generation  (SHG)  Wavelengths. 

(1)  Type  I  Phase  Match  Region.  The  plots  of  d^^  for  1.064/1.064/0.532  um  SHG  are 
shown  in  Figure  E-1  (KTP),  Figure  E-10  (RTP),  Figure  E-19  (KTA),  Figure  E-27  (RTA),  and 
Figure  E-36  (CTA).  The  d^g  plots  for  1.35/1,35/0.675  pm  SHG  are  shown  in  Figure  E-3  (KTP), 
Figure  E-12  (RTP),  Figure  E-20  (KTA),  Figure  E-28  (RTA),  and  Figure  E-36  ((H'A).  In  all  cases 
the  d^ff  goes  to  zero  at  the  <j)  =  0°  and  90®  extrema.  The  maximum  is  typically  near,  but  not 
exactly,  <J>  =  45°  due  to  the  varying  contributions  of  d,s  and  For  all  isomorphs  the  ^eff  values 
were  small  (see  Table  1 1);  RTA  and  KTP  had  the  largest  values,  but  even  they  were  only  on  the 
order  of  0.5  pm/V. 

(2)  Type  II  Phase  Match  Region.  Plots  of  d^g  for  1.064/1.064/0.532  pm  are  presented 
in  Figure  E'2  (KTT)  and  Figure  E-1 1  (RTP).  Only  KTP  and  RTP  are  shown  because  none  of  the 
arsenates  are  capable  of  Type  II  phase  matching.  For  1.35/1.35/0.675  pm  SHG,  d^g  plots  are 
shown  in  Figure  E-4  (KTP),  Figure  E-13  (RTP),  Figure  E-21  (KTA),  Figure  E-29  (RTA),  and 
Figure  E-37  (CTA).  The  maximum  d,g  was  found  where  (j>  was  at  a  minimum  value  and  0  was 
at  a  maximum  value.  The  minimum  occurred  at  (j)  =  90°  (which  is  also  where  the  minimum 
value  of  0  existed).  Differences  in  d,g  among  the  isomorphs  was  minimal  for  both  SHG  cases. 

(3)  Final  Comment.  When  phase  matching  was  possible,  RTA  had  the  largest  d,g  for 
SHG  in  both  Type  I  and  Type  II  cases.  The  exception  was  Type  II  1.064  pm  SHG  where  KTP 
was  the  winner  because  RTA  did  not  phase  match.  Type  11  phase  matching  d,g  values  are  always 
larger  than  those  found  for  Type  I  phase  matching;  about  six  times  larger.  However,  the  change 
ranged  from  nearly  24  times  (RTP  1.064  pm  SHG)  to  about  6.8  times  (RTA  1.35  pm  SHG). 

ii.  Optical  Parametric  Oscillator  (OPO)  Wavelengths 

(1)  Type  I  Phase  Match  Region.  Plots  of  d,^for  4.043/1.444/1.064  pm  are  shown  in 
Figure  E-5  (KTP),  Figure  E-14  (RTP),  Figure  E-22  (KTA),  Figure  E-30  (RTA),  and  Figure  E-38 
(CTA).  Plots  of  d,gioi  3.17/1.60/1.064  pm  are  shown  in  Figure  E-7  (KTP),  Figure  E-16  (RTP), 
Figure  E-24  (KTA),  Figure  E-32  (RTA),  and  Figure  E-40  (CTA).  In  all  cases  the  d^g  goes  to  zero 
at  the  (})  =  0°  and  90°  extrema.  The  maximum  is  generally  near,  but  not  exactly  at,  (})  =  45° 
due  to  the  varying  contributions  of  di^  and  ^24-  "Fhe  isomorph  RTA  for  4.043/1.444/1.064  pm  and 
3. 17/1. 60/1. OM  pm  has  values  of  0  which  are  significantly  larger  than  the  other  isomorphs 
maximum  d^g.  Between  the  other  isomorphs  the  difference  in  d,g  values  is  minimal.  Except  for 
RTA  values  (where  d^g  values  were  over  1.46  pm/V  at  4.043  pm  and  0.95  pm/V  at  3.17  pm) 
the  isomorph’s  d^g  values  were  fairly  small  (see  Table  1 1),  only  on  the  order  of  about  0.4  pm/V. 

(2)  Type  II  Phase  Match  Region.  The  plots  for  d^g  at  4.043/1.444/1.064  pm  are 
displayed  in  Figure  E-6  (for  KTP),  Figure  E-15  (RTP),  Figure  E-23  (KTA),  Figure  E-31  (RTA), 
and  Figure  E-39  (CTA).  For  3.17/1.60/1.064  pm,  the  d,g  plots  are  shown  in  Figure  E-8  (KTP), 
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Figure  E-17  (RTF),  Figure  E-25  (KTA),  Figure  E-33  (RTA),  and  Figure  E-41  (CTA).  Maximum 
dfff  was  where  9  was  maximum  and  <])  was  minimum;  minimum  was  where  9  was  minimum 
and  ({)  was  maximum.  The  difference  in  d^g  between  the  various  isomorphs  was  not  significant. 

(3)  Type  III  Phase  Match  Region.  There  are  no  plots  of  d^ffStt  4.043/1.444/1.064  pm 
since  none  of  the  isomorphs  phase  matched  at  this  wavelength  and  phase  match  type.  For  plots 
of  d^ffZt  3.17/1.60/1.064  pm  refer  to  Figure  E-9  (KTP),  Figure  E-18  (RTP),  Figure  E-26  (KTA), 
and  Figure  E-34  (RTA).  CTA  does  not  Type  III  phase  match  at  this  wavelength  combination. 
Maximum  and  minimum  d^^  values  occurred  as  at  Type  II  phase  matching.  The  largest  d^g  values 
were  found  in  KTP  and  KTA  and  the  difference  in  their  values  was  very  minimal.  The  difference 
between  the  other  two  isomorphs  (RTP  and  RTA)  d^g  values  was  minor. 

(4)  Final  Comments.  The  crystal  RTA  has  the  largest  in  the  Type  I  phase  matching 
regions.  The  crystal  CTA  has  the  largest  d,g  in  the  Type  11  phase  matching  regions.  At  Type  III 
phase  matching  (which  only  occurred  at  3.17  pm  OPO)  the  crystal  KTP  has  the  largest  d,g.  Type 
II  and  Type  IE  phase  matching  d^g  values  are  always  larger  than  those  found  for  Type  I  phase 
matching;  about  six  times  larger.  However,  the  difference  ranged  from  23  times  (CTA  4.043  pm 
OPO)  to  as  little  as  1.6  times  (RTA  4.043  pm  OPO). 


Table  11.  Isomorph  Maximum  d^giot  SHG  and  OPO 


Maximum  (pmA/) 

CRYSTAL 

SHG 

OPO 

TYPEI 

•TYPE  II 

TYPEI 

TYPEU 

TYPEI 

TYPEU 

TYPE! 

TYPEU 

TYPEUI 

KTP 

0.420 

3.268 

0.366 

2.603 

0.437 

2.387 

0.398 

2.412 

3.4%: 

RTP 

0.126 

2.815 

0.322 

3225 

0.342 

^751 

0.431 

2.768 

2.696 

KTA 

0.231 

— 

0.163 

2.781 

0.407 

^210 

0.174 

2.320 

3.148 

RTA 

0.46a 

— 

a494 

3J357 

Z320 

^413 

2.341 

aA 

0.216 

— 

0.226 

2.471 

0.134 

0.176 

— 

1.064/t.06«0.S32  |un 

1.35/1  J5A)  675  (un 

4.043/1.444/1.064  fun 

3.17/160/1.064  |un 

•  Shadod  regions  indicate  the  largest  value  in  each  column 


d.  Walkoff. 

Appendix  F,  Isomorph  Walkoff  Angles,  shows  walkoff  angle  as  a  function  of  angle  (j)  for  the 
thrf  e-wave  interactions.  Since,  in  a  biaxial  crystal,  both  eigenpolarizations  can  have  walkoff,  we 
show  true  values  of  walkoff  for  all  three  wavelengths.  Of  course,  for  efficient  wavelength 
conversion,  it  is  most  desirable  to  have  all  wavelengths  with  walkoffs  as  similar  as  possible. 

Walkoff  angles  in  the  KTP  isomorphs  never  exceed  3®.  This  corresponds  to  only  1  mm  of 
walkoff  in  a  2-cm  long  crystal  so,  in  all  cases,  walkoff  is  small  in  the  KTP  isomorphs.  However, 
with  judicious  selection  of  phase  matching  type  and  crystal  orientation,  walkoff  can  often  be 
made  to  approach  zero. 
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i.  Partial  Phase  Match  Regions. 

Interestingly,  for  isomorphs  in  which  phase  matching  did  not  exist  over  the  entire  xy  plane 
(i.e.,  (])  =  0°  to  90°)  and  where  the  value  of  0  approached  90°,  the  walkoff  angles  of  the  three 
wavelengths  converged  to  near  zero  values.  Type  II  1.064  pm  SHG  in  KTP  (Figure  F-2),  for 
example,  shows  this  convergence  effect.  However,  1.35  pm  doubling  in  KTP  is  phase  matchable 
for  all  (|)  angles  and  the  result  is  that  relative  walkoff  of  the  beams  never  gets  lea:  than  -2°  (see 
Figure  F-4).  RTP  (see  Figure  F-11)  and  CTA  (see  Figure  F-34)  were  the  only  other  isomorphs 
which  had  walkoff  convergence  at  1.064  pm  SHG.  No  isomorph  showed  this  effect  at  1.35  pm. 

For  OPO  interactions,  walkoff  convergence  occurs  only  for  RTP  (see  Figure  F-18),  KTA  (see 
Figure  F-26),  and  RTA  (see  Figure  F-34)  for  the  3.17/1.60/1.064  pm  triplet.  There  are  no  cases 
of  convergence  for  the  1.444  pm  signal  wavelength  case;  neither  OPO  case  displays  convergence 
for  Type  n  phase  matching. 

ii.  Walkoff  Variations. 

Among  the  isomorphs,  walkoff  for  the  Type  I  phase  match  region  for  RTP  (except  for 
1.064/1.064/0.532  pm  SHG  which  is  unremarkable  with  respect  to  the  other  isomorphs)  shows 
the  widest  walkoff  variation  by  far  across  the  phase  match  region  (see  Figure  F-12,  Figure  F-14, 
and  Figure  F-16). 

iii.  Walkoff  Similarities. 

(1)  Types  I  and  II  phase  matching  at  4.043/1.444/1.064  pm  OPO  for  KTA  and  RTA 
is  interesting  from  the  perspective  that  the  plots  of  the  Xj  walkoff  angles  appear  similar  except 
that  the  maximum  walkoff  angles  are  located  at  opposite  ends  of  the  xy  plane  with  respect  to 
each  other.  The  plots  show  a  rather  rapid  increase  in  difference  in  walkoff  (i.e.,  from  from  no 
difference:  0°)  between  and  X3  at  one  end  of  the  xy  plane  followed  by  a  consistent  and  steady 
decrease  in  the  difference  in  the  walkoff  angle  between  the  two  wavelengths  until  0°  of  walkoff 
is  again  reached  at  the  other  end  of  the  xy  plane.  For  instance,  in  Type  II  phase  matching  (see 
Figure  F-23  for  KTA  and  Figure  F-31  for  RTA),  the  maximum  amount  of  walkoff  between 
and  Xj  for  KTA  is  located  in  the  first  half  of  the  xy  plane  ((j)  =  0°  to  35°)  whereas  in  RTA  it  is 
located  in  the  latter  half  ((j)  =  65°  to  90°).  In  the  other  isomorphs  the  difference  in  walkoff  angle 
between  X,  and  and  one  of  the  other  wavelengths  (^  or  X.3)  is  minimal  throughout  the  xy  plane. 
The  walkoff  angle  differences  noted  are  pronounced  compared  to  the  other  isomorphs. 

(2)  For  Type  I  phase  matching  for  KTA  and  RTA  at  3.17/1.60/1.064  pm  OPO,  a 
similar  walkoff  phenomenon  occurs  as  at  4.043/1.444/1.064  pm  OPO.  See  Figures  F-24  and  F-32. 
In  this  case,  the  walkoff  for  KTA  between  X,  and  is  maximum  at  ({>=0°  and  approaches  a 
minimum  at  (}>=90°.  In  RTA,  the  reverse  occurs  (minimum  at  0=0°  and  maximum  at  0=90°);  the 
walkoff  variation  in  KTA  is  somewhat  less  pronounced  than  in  RTA. 
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7.  CONCLUSIONS/OBSERVATIONS. 


The  choice  of  crystal  for  any  particular  application  depends  upon  a  number  of  things 
(application,  availability,  cost,  wavelength  at  which  it  will  be  used,  temperature  considerations, 
absorption,  transmission  range,  etc).  However,  this  study  has  a  purely  theoretical  focus  and 
maximized  performance  is  considered  based  upon  phase  matching,  d,g^  walkoff,  and  transmission. 
Table  12  synopsizes  the  results  by  making  recommendations  for  the  crystal  to  use  at  what  phase 
matching  type  and  at  what  three-wavelength  interaction. 

a.  Crystal  Characteristics. 

Clearly  no  one  single  crystal  does  it  all  and  does  it  well.  However,  the  following  comments 
are  valid  for  making  recommendations. 

i.  SHG. 

(1)  Nonlinear  Coefficient.  Type  II  SHG  has  values  which  are  six  times  higher  than 
those  of  Type  I.  The  crystal  RTA  has  the  highest  value  for  SHG  except  at  Type  II  1.064  pm 
SHG  (where  KTP  has  the  highest  value).  KTP  has  the  next  highest  Type  I  d,ff  while  RTF  has 
the  next  largest  Type  II  d,ff. 

(2)  Phase  Matching.  The  crystals  RTP  and  CTA  are  highly  biaxial  (i.e,  the  value  of 
0  changes  considerably  throughout  the  xy  plane)  in  nature  compared  to  the  other  isomorphs.  KTA 
is  the  most  uniaxial-like  (i.e.,  changes  in  0  are  minimal  throughout  the  xy  plane)  of  the  isomorphs 
at  all  SHG  wavelengths. 

ii.  OPO. 

(1)  Nonlinear  Coefficient.  At  Type  I  phase  matching,  the  crystal  KTA  has  the  highest 
values;  the  next  highest  levels  are  with  the  crystals  RTA  (4.043/1.444/1.064  pm)  and  RTP 

(3.17/1.60/1.064  pm).  For  Type  n  phase  matching  the  highest  d,ff  values  are  with  the  crystal 
CTA;  the  crystal  RTP  has  the  next  highest  values.  For  Type  HI  phase  matching,  KTP  has  the 
highest  value  of  the  crystal  KTA  has  the  next  highest  value.  Types  n  and  HI  phase  matching 
d,ff  values  ore  2  to  4  times  higher  than  those  of  Type  I. 

(2)  Phase  Matching.  KTA  is  the  most  uniaxial-like  (i.e.,  ~  of  the  isomorphs  at 

all  Type  I  OPO  wavelengths.  RTA  is  the  most  uniaxial-like  for  Type  n  OPO  wavelengths.  KTP 
was  the  most  uniaxial-like  of  the  isomorphs  at  Type  III  phase  matching;  however,  not  to  the  same 
degree  as  either  KTA  and  RTA. 

iii.  General  Considerations. 

(1)  The  larger  |n,  -  n^l  differences  in  CTA  and  RTP  resulted  in  larger  variations  of 
0  over  the  phase  match  range  of  0  for  each  phase  matching  wavelength  and  type.  Variation 
among  the  other  isomorphs  was  minimal.  However,  variation  in  0  for  KTA  was,  overall,  the  least 
of  all  the  isomorphs.  In  other  words,  KTA  was  the  most  uniaxial-like  isomorph. 

(2)  Maximum  does  not  occur  exactly  where  walkoff  between  the  3-wave 
interactions  is  minimized  at  any  wavelength  or  phase  match  type.  Maximizing  d,^  while 
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minimizing  walkoff  (i.e.,  maximizing  signal  out  of  the  crystal)  necessarily  represents  a  tradeoff 
between  the  two  (d,ff  and  walkoff  angles).  However,  in  those  crystals  where  a  limited  range  of 
(j)  phase  matching  v^ues  exists,  maximum  is  very  close  to  where  the  difference  in  walkoff 
between  the  three  waves  is  nonexistent  or  extremely  minimal.  Ihus,  this  represents  a  case  where 
tradeoffs  between  walkoff  and  maximum  are  not  necessary.  However,  this  works  only  for 
specific  wavelengths  and  phase  matching  types.  In  this  study  this  situation  exists  only  at  Type 
II  1.064/1.064/0.532  pm  SHG  (KTP  and  RTP),  Type  II  1.35/1.35/0.675  pm  SHG  (CTA),  and 
Type  ni  4.043/1.444/1.064  pm  OPO  (RTP,  KTA,  and  RTA). 

(3)  Maximum  for  Type  I  phase  matching  is  basically  in  the  central  area  of  the 
phase  match  region.  For  Type  II  and  III  phase  matching  the  maximum  is  at  the  smallest  angle 
^  (and  largest  value  of  9)  at  which  phase  matching  occurs. 

iv.  Crystal  Transparency. 

Choice  of  a  crystal  is  not  usually  independent  of  the  wavelength  at  which  it  will  be  used. 
Therefore,  it  is  important  to  know  the  light  transmission  through  the  crystal  at  the  chosen 
wavelength(s).  Figure  16  and 
Figure  17  [13,  14]  provide  an 
indication  of  the  amount  of  light 
transmission  at  a  given  wavelength 
for  each  of  the  KTP  isomorphs. 

The  crystals  used  were  uncoated 
samples.  But,  because  their 
thicknesses  are  unknown, 
correlating  the  percent  of 
transmission  values  indicated  in 
each  of  the  figures  is  not  possible. 

Therefore,  the  percent  of 
transmission  between  figures  (and 
curves)  should  be  taken  as  a 
relative  data  item  only.  With 
appropriate  antireflection  (AR) 
coatings,  the  percent  of  light 
transmission  of  each  of  the 
isomorphs  can  be  made  nearly 
perfect  for  well  grown  crystal  samples,  at  least  in  the  “flat”  portion  of  the  transmission  region. 
It  is  important  to  note  where  the  transparency  of  the  crystal  begins  to  fall.  This  is  where  the 
amount  of  absorption  becomes  a  factor  for  consideration  since  losses  due  to  absorption  are  not 
compensated  for  by  AR  coatings. 

The  extended  transparency  in  the  infrared  of  the  arsenate  crystals  (KTA,  RTA,  and  CTA)  is 
clearly  an  advantage  over  the  phosphate  crystals  (KTP  and  RTP)  for  midinfrared  purposes  [13]. 
Note:  the  isomorph  transmission  curves  in  Figure  16  and  Figure  17  are  digitized  representations 
of  the  original  transmission  curves  which  were  found  in  hardcopy  figures.  These  digitized  curves 
were  then  replotted  here.  The  transmission  curves  for  CTA,  KTA,  and  KTP  were  obtained  from 
Reference  [13]  and  the  curves  for  RTP  and  RTA  were  obtained  from  Reference  [14]. 


Wavelength  (|xm) 

Figure  16.  Phosphate  Isomorph  Transmission  Spectra 
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(1)  SHG.  Since  each  of 
the  isomorphs  has  relatively  flat 
transmission  curves  at  the  SHG 
wavelengths,  this  region  of  the 
isomorph  transmission  curves  is  not 
displayed  in  cither  of  the  two 
figures  above.  Transparency  at 
SHG  wavelengths  in  this  study  is 
not  a  useful  discriminator  as  to 
which  crystal  to  select  or  reject. 

(2)  OPO.  For  the 
following  discussion  see  Figure  16 
and  Figure  17  for  references  to 
isomorph  light  transmission  curves. 

The  wavelengths  of  1.444  and  1.60 
pm  are  u.sed  with  1.064  pm  to 
generate  the  mid-infrared  (mid-IR) 
wavelengths  of  4.043  and  3.17  pm,  respectively. 

The  wavelength  of  3.17  pm  is  in  the  relatively  flat  transmission  curve  region  for  all  of  the 
isomorphs  (except  for  KTP)  which  indicates  that  the  arsenates  and  RTP  arc  good  candidates  for 
use  at  this  mid-IR  wavelength.  Beyond  this  wavelength,  however,  only  the  arsenates  are  good 
mid-IR  candidates.  At  3,17  pm,  the  KTP  crystal  indicates  about  18%  loss  and  for  RTP,  the  loss 
is  about  4%.  For  the  arsenates,  both  CTTA  and  KTA  experience  about  an  8%  loss  while  the  loss 
for  RTA  appears  minimal  at  best. 

At  4,043  pm,  the  two  phosphates  (KTP  and  RTP)  are  in  steep  decline  with  respect  to  their 
ability  to  transmit  light  (see  Figure  16),  So,  this  essentially  represents  the  maximum  mid-IR 
wavelength  at  which  they  can  (or  should)  be  used  (if  the  amount  of  absorption  can  be  tolerated). 
On  the  other  hand,  each  of  the  arsenates  still  has  a  nearly  flat  transmission  curve  at  this  point. 
Phosphate  losses  at  4.043  pm  arc  over  30%  for  KTP  and  around  17%  for  RTP— fairly  serious 
losses.  For  the  arsenates,  losses  OiC  around  8%  for  CTA  and  KTA  and  minimal  for  RTA. 

The  transmission  curve  for  KTP  shows  that  it  is  most  applicable  to  wavelengths  less  than 
2.75  pm  (the  flat  region  of  its  transmission  curve)  and  perhaps  a  small  region  around  2.83  to 
2.9  pm.  The  crystal  RTP  appears  applicable  to  those  wavelengths  from  about  3,2  pm  and  less  and 
perhaps  a  small  region  around  3.4  to  3.8  pm.  The  arsenates,  on  the  other  hand  appear  better 
suited  for  wavelengths  beyond  where  the  phosphates  arc  capable.  Both  CTA  and  KTA  appear 
well  suited  to  wavelengths  to  around  3.65  pm  and  RTA  to  around  4,6  pm.  A  5%  loss  from  the 
isomorph’s  “flat”  transmission  curve  region  was  used  as  an  arbiffary  demarcation  point  beyond 
which  absorption  losses  become  a  concern. 


Wavelength  (jim) 

Figure  17.  Arsenate  Isomorph  Transmission  Spectra 


b.  Conclusions, 

Table  !2  below  synopsizes  the  information  available  on  the  isomorphs  for  use  at  a  given 
wavclenr'.h.  pha.se  match  type  and  transmission  loss.  Since  there  is  little  utility  in  using  any  of 


34 


the  crystals  in  Type  I  phase  matching  only  the  isomorphs  with  die  maximum  are  provided  in 
the  table.  As  can  be  seen,  the  isomorphs  all  appear  to  perform  in  a  manner  similar  to  KTP. 

i.  SHG.  All  wavelength  selections  for  the  SHG  cases  were  in  the  spectral  regions  where 
losses  due  to  absorption  were  not  an  issue. 

(1)  1.064/1.064/0.532  pm,  Type  II.  Only  KTP  and  RTP  phase  matched  here.  Both 
perform  about  the  same.  However,  the  for  KTP  is  better  and  the  walkoff  is  about  three  times 
less  than  that  of  RTP.  Therefore,  KTP  is  the  best  option. 

(2)  1.35/1.35/0.675  pm,  Type  II.  All  isomorphs  phase  matched  and  all  performed  well 
in  this  region.  However,  RTA  was  slightly  better  in  terms  of  d,g  value  and  walkoff  (which  is  just 
over  half  that  of  the  other  isomorphs  except  for  CTA).  Although  CTA  had  by  far  the  smallest 
walkoff  (about  four  times  less  than  RTA),  it  also  had  the  smallest  value  for  d^g. 

ii.  OPO.  In  this  type  of  operation,  isomorph  transmissivity  becomes  a  major  consideration 
in  the  preferred  choice  of  crystal. 

(1)  4.043/1.444/1.064  pm.  Type  II.  Both  CTA  and  RTA  perform  very  well  in  this 
region.  The  for  CTA  is  slightly  higher  than  for  RTA,  its  walkoff  is  about  less  than  for 
RTA,  but  CTA  appears  to  have  higher  absorption  than  RTA.  The  transmission  characteristics  for 
RTA  show  it  to  have  a  relatively  flat  transmission  curve  out  to  about  4.25  pm  while  for  CTA, 
the  flat  portion  of  the  transmission  curve  extends  to  only  about  3.7  pm.  The  better  transparency 
for  RTA  at  4.043  pm  may  compensate  for  its  slightly  lower  with  respect  to  CTA.  Therefore, 
RTA  is  the  preferred  crystal  in  this  group  with  CTA  a  very  close  second. 

Although  the  phosphates  (KTP  and  RTP)  have  slightly  smaller  values  of  d^p  the  walkoff  for 
each  was  around  one-third  larger  than  that  of  CTA  (and  about  25%  larger  than  for  RTA).  Both 
phosphates  suffer  substantial  transmission  losses  which  makes  them  poor  OPO  candidates  for 
generating  4.043  pm.  The  crystal  KTA  has  the  smallest  value  of  d^g,  a  slightly  higher  walkoff 
than  RTA,  and  about  the  same  amount  of  transmission  loss  as  CTA.  In  short,  KTA  is  a  better 
choice  than  the  phosphates  but  is  not  quite  as  good  a  choice  as  either  CTA  or  RTA. 

(2)  3.17/1.60/1.064  pm,  Type  II.  The  isomorph  CTA  is  clearly  the  crystal  of  choice 
in  this  region.  It  has  by  far  the  smallest  walkoff  of  any  other  isomorph  and  transmission  losses 
are  negligible.  The  phosphates  have  decent  values  but  each  has  large  values  of  walkoff  (nearly 
3°).  KTP  has  significant  transmission  loss  which  makes  it  poor  a  OPO  candidate  for  3.17  pm. 
RTP,  on  the  other  hand,  has  a  small  amount  of  transmission  loss.  It  might  make  a  decent 
candidate  for  OPO  use  here.  RTA  has  a  slightly  less  d^g  value,  about  0.6°  more  walkoff  than 
CTA  and  negligible  transmission  loss;  KTA  has  the  smallest  d^^,  about  0.9°  more  walkoff  than 
CTA  and  about  the  same  transmission  loss.  In  short,  the  arsenates  are  the  better  choices  for  3.17 
pm.  Type  II  OPO  work  tlian  the  phosphates. 

(3)  3.17/1.60/1.064  pm.  Type  III.  The  arsenate  RTA,  while  it  does  have  the  smallest 
d^ff  value,  has  very  minimal  walkoff  and  negligible  transmission  loss.  KTP  has  much  to  much 
transmission  loss  to  be  a  viable  candidate,  even  though  it  does  have  the  largest  value  of  d^^.  The 
isomorph  KTA  has  small  walkoff  but  its  transmission  loss  cannot  be  ignored.  RTP  has  a  larger 
d^g  than  RTA,  four  times  the  walkoff  (though  it  is  still  fairly  small).  However,  it  has  a  small 
amount  of  transmission  loss  that  RTA  doesn’t  have.  Therefore,  RTA  is  the  isomorph  to  use  here. 
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iii.  General. 

The  arsenates  characteristically  have  a  smaller  walkoff  than  the  phosphates,  typically 
somewhere  around  0.8°  to  1.5°.  However,  this  represents  about  0.5  mm  in  a  2-cm-long  crystal. 
With  careful  consideration  and  design  (with  respect  to  crystal  size  and/or  pump  beam  size,  for 
instance)  walkoff  effects  can  be  minimized.  Where  KTP  does  not  have  the  largest  the 
difference  from  the  isomorph  which  does  is  typically  no  more  than  about  0.7  pm/V. 

The  bottom  line  is  that,  stricdy  from  a  performance  characteristics  basis,  as  long  as  the 
wavelengths  are  less  than  about  2.75  pm  there  appears  to  be  no  real  reason  to  switch  from  KTP 
to  another  isomorph.  At  the  OPO  wavelengths,  however,  RTP  and  RTA  are  good  choices 
primarily  because  their  losses  are  less  than  that  of  CTA  (RTP  has  about  half  of  the  loss 
experienced  by  CTA  while  RTA  has  negligible  loss).  Beyond  3.17  pm  only  the  arsenates  are 
good  choices  because  of  the  significant  transmission  losses  experienced  by  the  phosphates. 

Developing  the  necessary  expertise  to  grow  one  or  more  of  these  crystals  for  use  over  the 
choice  of  KTP  is  not  justified  given  the  indicated  performance  characteristics,  the  cost,  and  the 
necessary  time  provided  the  wavelengths  being  considered  are  less  than  2.75  pm.  Beyond  this 
is  a  region  in  which  KTP  experiences  severe  losses  and  cannot  adequately  compete. 

Funher  study  is  necessary  to  more  fully  characterize  the  isomorphs  RTP,  KTA,  RTA,  and 
CTA  for  applications  at  wavelengths  beyond  2.75  pm.  Other  differences  in  the  isomorphs,  such 
as  the  power  handling  ability  (where  bulk  and/or  surface  damage  occurs),  may  exist  but  are  not 
idendfied  in  (and  are  outside  the  scope  of)  this  study.  These  differences  could  prove  important. 
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Table  12.  Isomorph  Selection  Synopsis 
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- 

u 

HIH 

CTA 

3.145 

Bi 

29.85 

O' 

_ 

S 

RTP 

2.768 

WM 

50.74 

■■ 

4 

1 

KTP 

2.412 

48.56 

18 

1 

RTA 

2.341 

40.91 

3.17/1.60/1.064 

1 

i^n 

KTA 

2.220 

mm 

45.10 

mmm 

1 

III 

KTP 

3.496 

0 

O' 

18 

1 

KTA 

3.148 

BI 

4.48 

4.40 

- 

1 

RTP 

2.696 

11.33 

11.60 

4 

1 

RTA 

2.341 

mm 

0’ 

3.21 

- 

Notes: 

1.  Isomorph  selections  are  provided  according  to  decreasing  value  of 

2.  Because  the  Importance  of  Type  i  phase  matching  Isn't  signillcant  compared  to  Types  II  and  III,  only  the 
isomorphs  with  the  largest  devalue  are  listed. 

3.  Isomorphs  In  Types  II  and  III  which  are  predicted  to  be  the  best  performer  are  highlighted. 

Indicators: 

'  Wavelength  f-  >  vhich  the  walkoff  values  of  the  other  two  wavelengths  are  measured. 

-  Transmis^nit  ,04s  Is  minimal. 

♦  Preferred  selection  In  category 
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appendix  a 


Phase  Match  Calculations  Worksheet 
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CALCULATE  PHASE  MATCHING  ANGLES  WofkdieetOevelopedby  CiptOaleLFeaimoR 

for  Positive  or  Negative  Nonlinear  BiaxialAJniaxial  Crystals  wl/elos.  Writbt-PMtenan  afb  oh 

SellmdcfViiutf  CRYSTAL; _ 

S  :  =READPRN(seJlmicr  rxa)  - 

(222681  0.99616  0.21423  0.01369  \ 

1.97756  1  25726  0  20448  0  00865 

228779  1  20629  0.23484  0.01583  / 


ORIGIN  =0 

TOL  =0  000001 

t- — —  <»=  Degrees  to  Radiant  CoDversioa 

180 

(A-2a.b) 

WAVELENGTHS  SELECTED 

User  Provided  m  => 

X  ojHI  064 

<==  Enter  user  provided  values  (#I  and  112)  such  that  Xa  S  Xp 

(A-3a) 

User  Provided  lf>  => 

XpSI.60 

{A-3b) 

Xeaic  is  positive  valued 

^calc‘~(^a  “^(3 

when  Xot  <  Xp  => 

j  X^^  =3.1761194 

(A-3c) 

Assign  Wavelengths  to  satisfy:  Xl  2:  X2  >  X3. 

^  1  -"^calc 

^2  ”^P 

X  j.=X^ 

(A-3d.  c,  0 

X  ,  =3.1761 

Xj  -1.6000 

Xj  =  1.0640 

SET  Ihe  Phuc  Match  Category  with  which  you  desire  to  work  according  to  the  foltowlDg  CATECiORIES 
(NOTE:  SELECTING  A  PHASE  MATCHINGIS  NOT  REQUIRED!): 

1  -  SHG/SfM  2  -  DFM-32/OPA-32 

3  -  OFM-31/OPA-31  4  -  OPO 

cat:=i 

Sdlmeler  Equations :' 


Indices  of  Refraction 

ns(X,)  na^L  3) 

n.=  ny^X,)  07(^3) 

fU^X  jy  lu^X  2^  nr^X  3^ 

WRriEPR.N  ^indiccsd  =ni 


11.7577  1.7906  1  8027  \ 

1  7756  1.7982  1  8092 

1  8278  1  8655  1  8809  / 


(A4) 


(A-Sa,  b,  c) 


(A-5d) 
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SETUP  Indicts  of  Rtflractloa  for  Ftiase  Match  Equation  Use 


•2:=(»o.i)'^ 

•3  •  =  (“0,2)*^ 

(A-6a,  b,  c) 

bi:=(“i.o)‘^ 

b2:  =  (>>i,i)'^ 

b3:  =  (n,_2)'^ 

(A.«,e,0 

*  —  /  \*2 
'i-("2,oy 

=  2-  =  ("2.i)'^ 

*^3  ”("2,2)  ^ 

(A-6g,  h,  i) 

The  following  equiuoos  were  developed  using  the  1987  vticle  printed  in  Applied  Physics,  1987,  wntlen  by  Fshlen  and  Yao. 
SETUP  for  Phase  Matching  Calculations 

lot(0.^)  •=sai(e)-cos(o)  lcy(0.0)  :=sin(0)-sin(o)  kz(6.o)  :=co$(0) 

Setup  Each  Wavelength  for  Phase  Matching  Calculations 

2  2  2 

B  i(0,p)  :=-kx(0,<i)  •(bi-t-ci)-ky(0,(j)  (a  j^c  j) -kz(0,«)  -(a  j-t-b  i) 

2  2  2 

c  1(0.0)  :=la  (0.0)  b  I'C  i-t-ky(0,o)  a  j-c  n-lcz(0,o)  -ai-bj 

2  2  2 

B  2(0.0)  •=-kx(0,o)  ■(b2+C2) -''y(0.O)  ■(a2-t-C2)-''^(0''>)  •(•2+*>2) 

2  2  2 

C  2(0.0)  •=!« (0.0)  ■b2-C2  +  'O'(0.O)  •a2  C2-l-kz(0,o)  •a2'b2 

2  2  2 

B  3(0.0)  :=-lcx(0.O)  ■(b3-t-C3)-ley(0,o)  ■(a3-t-c  3) -kz(0,(>)  •(03-^03) 

'2  2  2 

c  3^0.0) -^la (0.0)  ■b3  C3-t-ky(0,o)  a3  C3-t-lcz(0,o)  ■a3  b3 


(A-7a.b.c) 


(A-8a) 

<.A-8b) 


(A-8c) 

(A-8d) 

(A-8e) 

(A-80 


Indices  of  RelVactlon  for  Phase  Matching  Angles  (Vs  1, 2,  &  3)  Column  I.  Slow  Ray,  Column2:  FastRay 


Wavelength  Dependent  Indices  Assignments 


-U.M 

(A-10a,b.c) 

>•1 

^2 

(A-lOd,  e,  0 

1 

^3 
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V 


AMipimcMU  for  Uie  A^rropriate  EqiMtioiw  to  Utc  in  SOLVING  fortkc  Phatt  Match ConWIoaa  forth* CATEGORY  Ihalwit 

An(0.*)  =if(CAT-1.n;^j<e.*),if(CAT-2,n;i1s(0.*).if(CAT-3.n;^,n;^j(e.0)))) 

=‘f(CAT-1.n;j^e.t),if(CAT-2.n;„(0.t).if(CAT-3.n;^{0.*).nx3j(0.*)))) 

ATin(0-*)  =if(CAT-1.n;^^0.t)..f(CAT-2.nxis{0.t).if(CAT-3.n;^2^0.*).n;^^0.*)))) 


(A-IU) 

(A-llb) 

(A-llc) 


SHGl(0.t)  =n;^,s(e.»)-n;^j(0.t)  SHG2(0.t)  =  n  ;^,,(0.*)  SHG3{0,0)  xis(0.0)  +  n;^2f<®**) 

(A-12a,b,c) 

OPOI(0.t)  =SHGI(0.t)  OPO2(0.t)  =SHG2(0.t)  OPO3{0,0)  =SHG3(0,0) 

(A-12dei) 

DFM32.1(0.t)  =njj^0.t)- n^2s(®-*)  DFM32_2(0.t)  =n;^3j(0.*)- DFM32_3(0»  =n;^3l(0,*)- n^2f(8.*) 

(A-t2&h.i} 

DFM31J(0.t)  =n;jf(0,t)- n^ls(0,»)  DFM31.2(0.t)  =n;jj(0.*)-n;i,(0,0)  DFM31_3(0,4)  =n;^|(0.0) -  n xi,(M) 

Ai(0,t)  =if(  CAT-1. SHGl(0,»),if(CAT«2.DFM32_l(0.t),if(CAT-3,DFM31_l(0,0),OPOl{0.*)))) 

Aii(9»  Mf{CAT-1.SHG2(0.t).if(CAT-2.DFM32_2(0.t),if(CAT-3.DFM31_2(0,*),OPO2(0,4)))) 

Aiii(0.*)  =if(CAT-1.SHG3(0,*).if(  CAT-2.  DFM32_3(0.t),if(CAT-3.DFM31_3(0.+),OPO3{0.*)))) 


(A-12jW) 

(A-12in) 

(A-12II) 

(A-120) 


SETUP  SOLV-E  BLOCKS  FOR  TYPE  I,  II.  AND  III  (m  •ppropriit*)  PHASE  MATCHING 
CAT  =  1.0000 


CRYSTAL; 


Solve  Block  Tor  Type  I  Phase  Matcl^C 

Given 

(A-I3) 

Ai(0.*)-Ati(9.*) 

Waveleneth  Denendent  Solution 

(A-13a) 

0<-  0>O  0>O 

2  2 

Constrain  to  Ouadnnt  1 

(A-I3b) 

•nutch(0-*)  =™D(0) 

Detennine  Solution 

(A>I3c) 

Solve  Block  for  Type  II  Phase  Matching 

Given 

(A-M) 

A  ii(0.t)-A  i-jj(0.t) 

Waveleneth  Denendent  Solution 

(A-I4a) 

0<-  0>O  ♦<-  t^O 

2  2 

Constrain  to  Quadrant  1 

(A-14b) 

»nutch(9-^)  =FWD(9) 

Detennine  Solution 

(A-I4C) 

Solve  Block  for  Type  III  Phase  .Matchin{ 

Given 

(A-I5) 

Waveleneth  Etenendent  Solution 

(A-ISa) 

0<^  020  ♦<-  #20 

2  2 

Constrain  to  Quadrant  1 

(A-I5b) 

=FIND(0) 

Determine  Solution 

(A-ISc) 
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Aj  =ATl(0i-^*)  Ago  =  1-6942  C  =A-n,(6i.#)  CgQ  =  1.6942  Ej  =A-nn(ei.#)  Eg^  =  1.6942 

(A.18..b.c) 


Phase  Matching:  Index  of  Figure  A-la  Figure  A-lb  Figure  A>Ic 

Refraction  of  Generated  Type  I  Type  II  Type  III 

vs  Applied  Wavelengths 

Vary  the  value  of  "4x_s<arf  below  (where  “X"  =1,2,  or  3  for  the  Xj  =3.1761  X2  =1.6000  X  3  =1.0640 

each  of  the  three  phase  matching  types)  to  locate  the  beginning  of 

the  phase  match  region  shown  to  end  at  ♦  =  90'  above  ^  .0000 


Type  I  Type  U  Type  111 

41_stait.=0  ♦2_sl*rt  =0  ♦3_start  =59  <===  Assign  start  of  phase  match  region  (A>19a,b,c) 

41  =4l_stait-5  42  =42_start-^  43  :  43_stait-^  «ndcSPyettt^dians  (A-19d,e.f) 

i_!!tart  =0  i_end  =90  i  =  i_stait..  i_eod  <===  Assign  test  ranee  for  fiaures  below  (A>20a,b,c) 


<«==  Assia 

n  ranee  as  6 

(A-21) 

=  ATi(0i.4l)  =1-7004 

Cj  =A-ni(0i.42) 

C._e^  =1.7004 

Ej  =A-nii(8i.*3) 

Ej  ^  =  1.695852 

(A-22a,b,c) 

=  Bi  end '^-7^1438 

D;  =A„(9i.42) 

D;e^  =  1.7250 

Ej  =Ai,i(ei,43) 

Fj  ^=1.695859 

(A.22d.e.f) 


Phase  Matching:  Index  of  Figure  .A-2a  Figure  A-2b  Figure  A-2c 

Refraction  of  Generated  Type  I  Type  II  Type  Ill 

vs  Applied  Wavelengths 
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Setup  for  D«te  FUt  for  Phase  Match  Anglei 


N  :0..90 

(A-27) 

'"'n 

phi^. 

theu,^  ='N-3 

(A-28a,  b,  c) 

theta  =06ig'9 

theta  [I  :  if/theta  j|  i;O,0ft«,j'9.O) 

N  \  N  ‘  / 

(A-28d,e.O 

theta  =006jsj9 

theu,[j^  =if(theu,ii^>O,0e0N-8.Oj 

(A-28g.h.i) 

"Matrix"  of  phase  match  region  data  values 

data  :augment(augment^augment(^augment(augmei«(augment(w  j, w 2), w 3), phi), theta  jJ, theta uj, theta mj 


DATA  FILE  FORMAT  {by  Column):  il.  12. 13.  ♦;  6  (I>:9  (II):  •  (tU) 
Saves  the  Phase  Match  Angle  pairs  for  Type  I,  II.  &  III  phase  matching 

WRITEPRN(PM.DATArta)  =daUi 


(A-29b) 


SCRATCH  PAD  AREA  .. 
43_stait' 59.0000 

9‘  =mnmch(5&V58.630'^) 
«‘=»'match<5&^'58-6299-^) 


0t -89.99997094 ‘deg 
et  =89.99997921  "deg 
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APPENDIX  B  Effective  Non-Linear  Coefficient  (d,^)  and  Walkoff  Calculations  Worksheet 
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CALCULATE  BIAXIAUUNIAXIAL  CRYSTAL  NONLINEAR  COEFFICIENT  (d,„) 
and  WALKOFF  (p)  FOR  POSITIVE/NEGATIVE  CRYSTALS 


OBTAIN  DATA  FROM  DATA  FILES 

Murix  :=REA0PRN(PM_DATA  rxa)  <1  ijk  ■=READPRN(dijk  n  .=READPRN(mdicc*d  lu)  (B-la,  b.  c) 


IDENTIFY  ORIGIN  FOR  WORKSHEET  MATRICES 
ORIGIN  :=o 

d|j|(  for  the  crystal 


0.00 

0.00 

0.00 

0.00 

1.90 

0.00 

0.00 

0.00 

0.00 

3.60 

000 

0.00 

2.30 

3.80 

15  80 

0.00 

000 

0.00 

(B-2) 

Indicts  of  Rtfractlon  (X  by  column  •  Itft  lo  right, 
axts  by  row  -  x,  y,  and  z) 

h. 757675  1.790561  1.802654  \ 

n=  1.775611  1.798208  1.809211 

\  1.827846  1.865507  1.880939/ 


Select  Positive  or  Negative  Crystal  Type  with  a  positive  or  negative  Number.  cr  :=- 1  (B-s) 

NOTE;  Unless  properly  selected,  the  calculated  value  for  and  waikoll  angle  may  not  be  correct. 

PHASE  MATCH  angles  data  flit  LEGEND  [Matrix],  (by  column):  X2,  &  X3;  4>;  6(1);  6  (II);  6  (III) 

SELECT  AN  ANGLE  6  FOR  WHICH  AN  ANGI.F  fl  EXISTS  (lor  Tvoe  I.  II.  or  III  Selected  Belowl 


i:=o 

(B-4) 

X  (Matrix'^ '  ^)j  (Malrix'^^^))  j 

(B-5a) 

^l’"^0.0  ^2”^0.1  ^3"^0.2 

Xj=3.176  X2  =  1-6  X3  =  1.064 

(B-5b.  c.  d) 

Selects  the  Phase  Matching  Angle  Pair  from  the  data  file 

TYPE  1  TYPE  II 

TYPE  111 

I  ;=augineal(M»irix'‘^^  .Matrix*^^)  n  ■.=augmeni(Matrix*^^^  .Matrix''^'*) 

n  ■.=augmeni(Matra^^^  .Matrix*'^^) 

(B-6a.  b,  c) 

kfantlfy  the  Tvne  of  Phase  Matching  Desired  ( f  s  Type  1, 2  =  Typ»  II,  3  s  Typ»  f/i) 

PMId:=1 

Type :  =if(  PMIDSl .  I.  if(  PMIDS»2 .  n .  in ) ) 

Ai:=(Type‘'°^),deg  9:=Ai  0  =40.3468*deg 

Pj  ■=(Type‘^*^)i-deg  (!i:=Pj  d  =  0*dcg 

SELECT  Phase  Match  Category  for  Desired  Nonlinear  Coefflclent(di.Hl 
1  -  SHG/SFM  ("Xs”)  2  -  DFM-32/OPA-32  ("Xi") 

3  -  DFM-31/OPA-31  {"X2")  4  -  OPO  ("Xl"> 

CATEGORY  =4  (B-8) 


(B-7a) 

(B-7b) 

(B-7c.d) 

(B-7o,  0 
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Th/s  r0glon  b*glnB  th»  ct/cu/stfon  lor  tho  oHaeUvo  non~llnoar  eootltdont  (dtu). 
Thoao  coleulaUons  an  NOT  llmitad  to  tha  prinelpla-axia. 


This  is  th«  solution  lor  tha  OPTIC  AXIS.  Q 


n.= 


“2.0 

("1.0)^  “  (“0.0)^ 

"1.0 

.("2.0)'"  ("0.0) 

*^.1 

("1.1)  ^“(“0.1) 

"1.1 

.("2.1)  “("0.1) 

"2.2 

("1.2)^- ("0.2)^ 

“1.2 

.("2.2)^“  (“0,2) 

■i.s' 

'.s' 

-1.5' 

. 

This  angle.  O.  is  measured  Irom  the  s-axis  to  the  optic  axis. 
NOTE,  the  10  ‘18  value  is  intended  to  prevent  tan2S  from 
experiencing  a  singularity 


131.107  \ 
19.167  -deg 
17.333  / 


(B-9) 


This  Is  used  to  calculate  the  polarization  angle,  “fi”.  This  angle  is  formed  by  the  direction  of  the  slow  ray  electric  displacement 
unit  vector,  e^,  from  the  surface  of  tlte  plane  formed  by  the  wavefront  vector  (A)  and  the  z-exis 


cot(6)-fln(2o) 


“‘(^o.o) 

2  2  2  2  2 
■sm(0)  -coi(0)  'COt(0/  d-nc(o) 

coi(0)  •nn(2-o) 

c<x(Oi_o] 

510(0)  —  c<>s(0)  ■co5(<l)  d-sm(o) 
cos(p) '501(2  0) 

sm(6)^  — c£k(0)  'CO((d}^4-dai,p) 


5  njj  l-if(co»(0)-«in(2'p)tSO.O,if(co«(0)-an(2'0)>O,O,tt)) 


5-  = 


(suii(un26o_o).Hfivl^‘JS 
(sun(t*n25i  o).f.6^j)  o.5 
(•trwi(ua252  o) -t-5^j)  0.5 


/O 

ui]2S2  I  0 

lo 


5|idj=.0 


cos(6).na(2-^)  >0 


(0.000  \ 
0.000  rdej 
0.000/ 


(B-lOa) 


(B-10b) 


(B'lOc) 


The  following  are  the  solutions  for  the  SLOW  (b1 )  rays  and  the  FAST  (b2)  rays.  These  solutions  represent  the  UNIT 
VECTORS  of  the  electric  displacement  vectors.  Columns  represent  Xi,  Xj,  &  X3,  respectively. 

b  ;=-coj(0)-un  W  -““(^o.  0)  -  “»(♦)  -“(fio.o) 
d  ■=-cot(0)  no($)  co»^8j  (})-cos(6)  im(8j  q) 
f :=-co»(e)-tio(^)  eos(62_Q)  - eos(0) -00^82  0} 


ii“-cos(8)-co$  ($)  •cos  (®o.o)  +  sin(^)-siny^  q) 
c  .=-cos(0)-co$  (0)  •cos  (8i,(j)+$in((>)-$m(8,.o) 
c  •^cos(0)-co«  (0)  ‘COS  (52_o)-t-an(«)-iin(62  0) 


as  •=-coi(o)  coi(o)-sui^8o  0)  -  «>“(p)  “*(fio.o) 

cc  cos(o) ■cos(o) -sm  ~  aa  ((j)-coi^S,  0) 

cc  •--cos{6)  005(0)  sui  (Sj  0)  -  «n(o)  c«(62.o) 

/a  c  e 

:|  b  d  f 

^sin(0)  imsICq  q',  sin\,0)  q)  sin(0)  coj(82  q) 


bb  =-cos(0)  tin(o)  $in(^SQ  q)  d-coa(<))-coa^8Q_Q) 

dd  ■&n  (0)  •sio  (8i_o)-f-cot(«)  cof(8i  0) 

ff '  =-«»(0) -sm  (<i)-iin(82  0)  -t-co*(«) •001(^82  0) 


b  t  .= 


b2'  = 


bb 


CC 

dd 


ee 

flf 


\sn(0)  ■sin  (®o.o)  sin(0)‘sm  (Sl.o)  00(0)  •sin  (52.0)1 


(8-1  laA) 
{B-1lc,cl) 
(B-II0.O 


(B-ilg.h) 

(8-l1i,i) 

(B-llk.O 


(B-lim.n) 
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Thl«  «>ctlon  c«ietil«f « th>  •l>cWefltld-n  two  aoliriiatlon  comaontnM.  Xi .  Xa,  and  X3  (top  to  bottom),  respectively. 
Column  1  is  (or  the  ‘SLOW  rays  (e^);  column  2  calculates  the  'FAST*  rays  (®2)- 


f  ^  ‘0.0)^ 

{‘’'1.0)^ 

(‘>•2.0)^' 

0.5 

'(‘‘^O.o)^ 

(‘>2..0)^ 

05’ 

0 

(“..0)' 

h.o)‘  . 

.  Ko)^ 

(“i.o)' 

("2.0)  . 

0.5 

0.5 

(‘'•o.t)' 

('’'i.t)' 

(^20.,)^ 

h2.y 

(“1,.)' 

(“2.1)'. 

.  ("0,.)' 

I**  ‘0.2) 

('»t.2)^ 

0.5 

'(•>20.2)^ 

(‘’2..2)^ 

(^^2,2)1 

0.5 

.[  (“0.2)^ 

("2.2)** . 

.  (“0.2)”* 

r  i 

(“1.2/ 

(“2.2)^  . 

This  section  calculates  the  UmT.VECTORS  (or  the  oolariratlon  components  for  the  electric  fields.  The  'SLOW*  rays  are 
calculated  by  "al the  ‘FAST*  rays  by  *32.*  The  rows  represent  the  x.  y.  and  2-axis  components,  respectively,  the  columns 
represent  Xi.  X2  & X3  (let!  to  nghl).  respectively. 


**  *0.0 

cr 

0 

^  *0.2 

^o.o('’o.o)^ 

p^i.oN.i)^ 

''2.o('’o.2)^ 

*>*1.0 

‘■>1.1 

‘**1.2 

^o.o(“i.o) 

^l.o("t.l)^ 

*^2,o("l.2)" 

>2.0 

‘■*2.1 

‘’*2.2 

.*'’0.0‘(“2.o)^ 

'^1.0-("2.lf 

2 

’''’2.0' ("2.2) 

'’2o.O 

cr 

0 

'’20.2 

‘’’’0.1  (“0.0)^ 

P‘’l.l("0.lf 

‘’*’2.1' (“0,2)^ 

'’2|.0 

'’21.1 

‘■21.2 

‘^o.i(“i.o)^ 

‘’'*2.l("l.2)^ 

*’22,0 

'’22.1 

'’22.2 

,^'o,r(“2.o)^ 

'^i.r(“2.i)^ 

'^2.^(»2.2)^ 

(B-13a.b) 


Categories:  CATEGORY  =4 

1  -  SHG«FM.  2  -  OPO-32/OFM-32 
3  -  OPO-31/OFM-31 . 4  -  OPO 


Crystal  Tvoe  (Positive  or  Negative 
Number  indicates  Type) 

CT=-1 


Phase  Match  Type 
PMlDal 


Identity  which  column  to  use  tor  the  specitic  phase  match  category  selected 

C  i  :=U(CATEGORY=1.2.if(CATEGORY=2,0.i((CATEGORYS3.l.2))) 

C,=2 

(ai) 

(B-I4a) 

C  j :  =if(  CATEGORY*! .  0,  if(  CATEGORY*2 . ! ,  if(  CATEGORY*3 .0.0))) 

Cj=o 

(aj) 

(B-I4b) 

C  k  : =a(  CATEGORY*! , ! .  a(  CATEGORY*2 , 2 ,  a(  CATEGORY*3 .2.1))) 

Ck  =  ! 

(aiO 

(B-14C) 

Setup  the  a,,  aj.  and  a,,  vectors 

Type  1.  It,  and  III 

*  j  ■.=if(rT>0.«2^^  '^) 

■=if(cT>0.»l^^  .•2'^^ 

Type  1 

Type  ll/lll 

»lj  =.f(cT>0.»l*^^*^.42'^^*^) 

aij  •-if(cT^0.a2  *  ,*!  *  )  »j  =if(PMlD>!  ,a2j,a  ij) 

(B-I5a.b) 


(8-1  Sc.  d.e) 
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SelUD  the  aja,,  column-vectors 

IXBftJ  IvatlMM 


*  >JO.O  *ko,0 

*2jo.O*ko.O 

*'ii.o*''i.o 

*2jl.0*kt.0 

*  ')2.0  “ki.o 

*2)2.0  *>'2.0 

*  ‘)l.o“k2,o'^®  02.0*kl.O 

»2jk  ■  = 

*2jl.0*k2.O^*2j2.0’*kl.0 

*  'JO.O  * >'2,0'^*  ')2.0'*k0.0 

*  2jo .  o'*  k2 . 0  ■*■  *  2j2  .  o'*  ko .  0 

.*  'io.o‘kt.0'''*  'J1.0*ko,0 

.*2jo,0'*kl.0‘*"*2jt,o'*ko.0. 

Sglact  the  eporopflate  a)ai>  column  vector 
•  jk  ■  =if(PMn3>  1 .  •  2jk .  •  1  ,it) 

Effective  Nonlinear  Coefficient  PMID  =  1 

eir  I*  i  (<*  ijk  «jk)|  <*  eff  =  2-33754626 

WALKQFF  Angle  tot  Calculatlona  (coiumni :  SLOW  rays,  column  2:  FAST  rays) 
By  row.  starling  at  the  top  is  the  calculation  lor  the  Xi,  X2.  &  2.3,  respectively 


p=Rd 


'o.o\ 


"0.0 


"1.0 


"2.0 


Pb, 


0.0 


'0.1 

"0.1 


bl 


1..1 


"l.t 


2  'h 

‘’'2.1' 


a 


"’l.O 

h.2\ 


"0.2 


■1.2 


"2.2 


Pb 


2.0 


^OA 

ri.ii 

\^aI 

1 

p^i.t 

^21.2! 

buii 

\“2.2i 

"IM 

Pb. 


2.1 


Select  correct  welkoff  angles  according  to  phase  match  caleoaiv 

v:=il2 

rt 

m  1  :=if^CATEGORY=2.if(PMlDSo.po_i.Po.o)-‘^(P‘'<®-°'Po.O'Po.l))''' 

w  2  :  =if(CATEGORYS:3 . if (PMED^O. P  j .  j . P|  .0)  •  if(PMID>0. p, . q. P 1 . 1))  v 
vw>  3  :  =if^PMtD2o. p2  (.P2  o)"'' 

PMID  =  1  ORIGIN  =0  CATEGORY  =4  CT  =“1 

Create  the  Matrix  for  Writing  a  Line  of  Data  to  the  Data  File 
d«u.=^X|  Xt  X3  8v  ov  djff  TO|  W02  *03^ 

Writes  Data  to  the  Data  File  Specified 
APPE.S'DPR.V^pmold  =d2Us 


(B-16C) 

(B17) 


(B-10) 


(B-19a) 

(B-19b) 

(B-19C) 

(B-19d) 


(B-20a) 


(B-20b) 
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CONSOUQATeP  iMEO.R!VIATION 


Etoctric  Displacement  Unit  Vectors,  b1  •  ‘slow*  ray,  b2  •  *tast‘  ray 


b  1  = 

'-0.76213977  -0.76213977  “0.76213977' 

0  0  0 

b2  = 

^0.64741252  0.64741252  0.64741252  i 

\ 

Polanzation  Unit  vectors,  at  >  'slow*  rays,  a2  <i  'last'  rays 


^.78640149  -0.78751432  -0.78841497\ 

jo  0  0\ 

al  = 

0  0  0 

tO.61771571  0.61629635  0.61514375 ; 

a2  = 

1  1  1 

,0  0  0/ 

Polarization  Components  for  tha 
Electric  Fields  (cols  =  "slow'*  & 
"fast"  respectively) 

(0.31369897  0.31717898  \ 
0.30185472  0.30925747 

029747807  0.30550741  / 


Polarization  Unit  vectors 

ai  >  'Generated  Wavelength*,  aj  and  ak  are  the  'Applied*  Wavelengths 


/-0.78841497\ 

l°\ 

1  h 

.i=  0 

•j=  Pj 

1  *jk  = 

\0.6t5t4375  / 

\oi 

\ol 

Column 

vector  Irom  'applied*  Wavelengths  =» 

Tensor 

/O  0  0  0  1.9  0 

4ijk=.  0  0  0  3.6  0  0 

\2.3  3.8  15.8  0  0  0 

Polarization  Vector 

4_ej»k:=dijk«jk 


d_sjak 


a  i-d_ii»k=  2.33754626 


Indices  of  Retraction  (A.  by  columns) 

11.7576751  1.7905609  1.8026541  \ 

1.7756109  1.7982079  1.8092106 

1.8278463  1.8655065  1.8809392/ 


WalkoK,  column  1  -  'slow*  rays 
column  2  -  'last*  rays.  X  by  rows) 


(2.19728498  0\ 
2J0062428  0  I ’deg 
2.38443356  0/ 


Xi  X2  X3 
P  r”(Po.o  Pi,o  P2.1)  '' 

P  n  -“(PO.!  Pl.O  P2.l)’'' 
P  in  -“(Po.o  Po.i  P2,i)‘'' 


WalKoH  Angle  in  degrees,  by  Type 
p  j=(  2.19728498  2.30062428  0  ) 

Pll={0  2.30062428  0  ) 
p  m  =(  2.19728498  0  0  ) 


Mtscelleneous  Information 

^  =  0‘deg  e  =  40.3468*deg 
X  j  =3.t761lX<2  =  1.6 
X3  =  1.064 

CATEGORY  =4 
1-  SHGTSFM,  2-OPA/DFM-32 
3-OPA/DFM-31.4-OPO 


Phase  Match 
Type: 

PMnJ  =  ' 

Angle  ol 
OPTIC  Axis 
from  z-axis 


Crystal  Type 
(+/-  number 
indicates  type) 

CT=-1 

/0.54291921  \ 
=  033452255 
\0.302S2184  / 


Ellectivs  Nonlinear  Coefficient 
d  eff  =  233754626 

Worksheet  Waltolf  values 
(according  to  phase  match 
type  selected) 

wo  1  =2.19728498 
wo  2  =2.30062428 
wo  3  =0 


NOTES: 

1 .  This  worksheet  treats  uniaxial  crystals  as  a  simple  case  ol  a  biaxial  crystal,  hence  3-axes  vs  2  (i.e..  ooe  vs  oe) 

2.  For  negative  crystals,  the  deliniton  ol  'slow'  and  'fast*  rays  are  'reversed*  from  what  they  mean  for  positive  crystals 

3.  Tha  elfactive  nonlinear  coeffioent.  cgh  is  same  lor  SHG/SFM  and  OPO  categories  (i  e .  categories  1  and  4).  Likewise,  the  two 
OPAfDFM  cases  (calegones  2  and  3)  have  the  same  cgh. 
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APPENDIX  C  KTP  Isomorph  Phase  Match  Ranges 
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List  of  Tables.  Appendix  C 


Table  C-1  Phase  Match  Range.  KTP .  54 

Table  C-2  Phase  Match  Range.  RTP .  55 

Table  C-3  Phase  Match  Range.  KTA  .  56 

Table  C-4  Phase  Match  Range.  RTA  .  57 

Table  C-5  Phase  Match  Range.  CTA  .  58 
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Table  C-l  Phase  Match  Range,  KTP 


-  KTP- 
Wavelength  (pm) 

Phase 

Match 

Type 

Phase  Match  Range 

8,(|)  (Degrees) 

B 

6 

4> 

A0 

A(t> 

1.064 

1.064 

0.532 

I 

49.41.  39.20 

0,90 

90 

II 

90.00,  68.67 

26.0,  90 

21.34 

64 

1.35 

1.35 

0.675 

I 

40.25,  29.57 

0.  90 

10.68 

90 

II 

49.11,  39.42 

0,90 

9.70 

90 

4.043 

1.444 

1.064 

I 

39.24,  29.06 

0.  90 

10.18 

90 

II 

44.33,  35.47 

0,90 

8.86 

90 

III 

;•  . 

• 

■1 

H 

3.176 

1.60 

1.064 

I 

37.62,  27.03 

0.  90 

10.58 

90 

II 

44.85,  35.96 

0,90 

8.89 

90 

III 

77.36,  65.90, 

0.  90 

11.46 

90 

5' 


Table  C-2  Phase  Match  Range,  RTP 


-  RTP  - 
Wavelength  (jim) 

Phase 

Match 

Type 

Phase  Match  Range 

9.(J)  (Degrees) 

mm 

8 

Ad) 

1.064 

1.064 

0.532 

■■ 

55.83,  32.17 

0,90 

23.66 

90 

II 

90,  64.62 

41.70,  90 

25.38 

48.3 

1.35 


1.35 


0.675 


I 

47.12,  18.82 

0. 

90 

28.31 

90 

II 

65.95,  43.91 

0. 

90 

22.05 

90 

1.064 

I 

47.61,  19.74 

0.90 

4.043 

1.444 

II 

52.79,  28.08 

0,90 

m 

27.88 

24.71 


II 


0,90 


0.  90 


43.35,  90 


30.08 


24.49 


23.32 


90 


90 


44.65 


3.176 


1.60 


1.064 


III 


45.85,  15.77 
53.16,  28.67 
90,  66.68 


Table  C-3  Phase  Match  Range,  KTA 


-  KTA  - 
Wavelength  (pm) 

.  ^l„.- 

JU 

mm 

1.064 

1.064 

0.532 

1.35 

1.35 

0.675 

4.043 

1.444 

1.064 

3.176 

1.60 

1.064 

Phase 

Match 

Type 


Phase  Match  Range 
0.(j)  (Degrees) 


51.47,  44.69 


37.30,  32.98 


45.72,  34.82 


A6  A6 


0,90  6.78  90 


I 

40.76.  33.79 

0,90 

6.97 

90 

11 

67.67,  54.67 

0,90 

7.00 

90 

0,  90  4.32  90 


0,  90  10.90  90 


[i] 


Table  C-4  Phase  Match  Range.  RTA 


APPENDIX  D 


KTP  Isomorph  Nonlinear  Characteristics 
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Table  D- 1  Nonlinear  Characterisiics,  KTP 
Table  D-2  Nonlinear  Characteristics.  RTF 
Table  D-3  Nonlinear  Characteristics.  KTA 
Table  D-4  Nonlinear  Characteristics.  RTA 
Table  D-5  NonUnear  Characteristics.  CTA 


Table  D-1  Nonlinear  Qiaracteristics.  KTP 


1  Starred  (*)  walk-ott  angle  indicates  the  walk-otf  angle  from  which  the  other  two  are  rel^renced 


2  All  measurerrients  made  where  d,„  was  at  maximjm  value 


3  SHG  Measurements  made  for'  1  064/1  064/0  532  and  1  35/1  35/0  675 


4  OPO  Measurements  made  lor  4  043'1  444/1  064  and  3  176/1  60/1  064 
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1  Starred  (*)  walK-oH  angle  indicates  the  walk-otf  angle  from  which  the  other  two  are  referenced 
2.  All  measurements  made  whe:e  d^  was  at  maximum  value 


3  SHG  Measurements  made  for-  1  064/1  064/0  532  and  1  35/1  35/0  675 


4  OPO  Measu'-emen's  made  for  4  043/1  444/1  064  and  3.176/1  60/1  064 
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Table  D-3  Nonlinear  Qiaractcristics.  KTA 


1  Starred  (*)  walk-ctf  angle  indicates  the  walk-otf  angle  from  which  the  other  two  are  referenced 

2  All  measurements  made  where  d.,  was  at  maximum  value 

3  SHG  Measurements  made  for  1  064/1  064/0  532  and  1  35/^  35/0  675 

4  OPO  Measurements  made  for  4  043/1  444/1  064  and  3  176/1  60/1  064 
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Table  D-4  Nonlinear  Characteristics.  RTA 


1  Starred  (*)  walk-ott  angle  indicates  ttie  waik-oH  angle  from  which  the  other  two  are  referenced 
2.  All  measurements  made  where  d.,  was  at  maximum  value 

3  SHG  Measurements  made  for  1  064/1  064/0  532  and  1  35/1.35/0  675 

4  OPO  Measurements  made  lor  4  043/1  444/1  064  and  3  176/1  60/1  064 
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Table  D-S  Nonlinear  Oiaracteristics.  CTA 


1  Starred  (')  waJk-ott  angle  indicates  the  waik-otl  angle  from  which  the  other  two  are  referenced 

2  All  measurements  mar'a  whore  d^  was  at  maximum  value 

3  SHG  Measurements  made  (or  1  064/1  064/t)  532  and  1  35/1  SS-X)  675 

4  OPO  Measurements  made  for  4  043/1  444/1  064  and  3  176/1  60/1  064 
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APPENDIX  E 


KTP  Isomorph  Nonlinear  Coefficient 
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Figure  E-20  versus  Phase  Match  Region,  1.35/1.35/0.675  pm.  Type  I,  KTA 
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Figure  F-5  Walkoff.  OPO,  4.043/1.444/1.064  lam.  Type  I.  KTP 


Figure  F-6  V^alkolf.  OPO.  4.043/1.444/1.064  urn.  Type  II.  KTP 


96 


Walkoff  (Deg) 


Walkoff  (Deg) 


10  20  30  40  50  60  70  80  00 

♦  (Deg) 

Figure  F'l2  Walkoff.  SHG,  1.35/1.35/0.675  pm.  Type  I.  RTF 
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Figure  F-13  W.ilkoff.  SHG.  1.35/1.35/0.675  pm.  Type  11.  RTF 
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Figure  F-14  Waikoff.  OPO,  4.043/1.444/1.064  pm.  Type  I.  RTP 
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Figure  F-15  Waikoff.  OPO.  4.043/1.444/1.064  urn.  Type  II.  RTP 
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Figure  F-23  Walkoff.  OPO.  4.043/1.444/1.064  pm.  Type  II.  KTA 
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Figure  F-26  Walkoff.  OPO.  3.17/1.60/1.064  pm.  Type  III,  KT.-\ 
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Figure  F-28  Walkoff.  SHG.  1.35/1.35/0.675  pm.  Type  1.  RTA 
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Figure  F-29  Walkoff.  SHG.  1.35/1.35/0.675  pm.  Type  II.  RTA 
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Figure  F-37  Walkoff.  SHG.  1.35/1.35/0.675  pm.  Type  II.  CTA 
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Figure  F-40  Walkoff,  OPO.  3.17/1.60/1.064  pm.  Type  I.  CT.a. 
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Figure  F-41  Walkoff.  OPO.  3.17/1.60/1.064  pm.  Type  II,  CTA 
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DEPARTMENT  OF  THE  AIR  FORCE 

WRIGHT  LABORATORY  (AFMC) 
WRIGHT-PATTERSON  AIR  FORCE  BASE.  OHIO 


21  March  1994 


MEMORANDUM  FOR  DTIC/OCC 

CAMERON  STATION  BLDG  5 
ALEXANDRIA  VA  22304-6145 


FROM:  WL/ELO  Bldg  22B 
2700  D  St  Ste  2 

Wright-Patterson  AFB  OH  45433-7405 


SUBJECT:  ChangesAJpdates  to  WL-TR-93-5093,  Three-Wave  Nonlinear  Interactions  in  KTi0P04  (KTP)  and 
KTO  Isomorphs  (AD  Number  A274079) 


As  with  many  publications,  there  are  inevitable  errors  in  the  final  product  which  need  correction-and  this  report 
is  no  exceptimi.  Following  are  changes  and  COTrections: 

Note  The  distinction  between  “ray”  and  “v/ave”  is  somewhat  loosely  treated  in  the  report. 

Properly  caisidered,  treat  as  a  reference  to  a  “wave”  and  S  as  a  reference  to  a 
“ray.”  Specifically,  k  refers  to  the  wave  normal  vector  (direction  of  propagation)  and 
S  refers  to  the  ray  vector  (direction  of  energy  flow). 

Page  4  Figure  3  incorrectly  shows  the  angle  ([>  as  being  referenced  from  the  y-axis.  Modify  the 

figure  to  indicate  reference  as  being  from  the  x-axis. 

Page  6  Figure  4  was  obtained  from  reference  [4].  Credit  for  this  figure  was  inadvertently 

removed  from  the  Final  Report 

Page  10  Last  paragraph.  Replace  the  sentence 

The  value  of  d,n  is  d^endent  upon  the  direction  of  the  1^...  we  will  use  the  electric 
displacement  vectors,  tr((oJ  to  calculate  aj,  a,,  and  a^). 

with 

The  value  of  d,„\s  dependent  upon  We  do,  however,  know  that’^JJ?  and  that"^!^.  We 
also  know  that  D  and*^  are  related  in  terms  of  the  polarizability  tensor.  Therefore,  we  can 
use  the  electric  displacement  vectors, "5, (coj,  to  calculate  aj,  a,,  and  a*  of  eqn  16. 

Page  12  Change  the  indices  of  refraction  subscripts  in  eqn  23  from  “1,2  and  2”  to  “x,  y,  and 

z”,  respectively. 

Page  14  2nd  Line.  Change  “...the  same  as  {iijoij}  to  “...the  same  as  {a,J(Og)} 

Paragraph  3.b.ii,  2nd  line.  Change  “...which  exactly..."  to  “...which  are  exactiy..." 

Page  15  Paragraph  3.b.vii.(l).  Change  “...exactly  the  value"  to  “...exactly  the  same  value" 


Page  19  Table  2,  last  cell.  The  Sellmeier  equation  should  read: 

n?  =  Ai  +  (Bi/[l  -(CA)"]-D.A,‘ 

Page  30  Table  11.  Orange  d^f  for  KTP,  0.675  pm-SHG,  Type  II  from  2.603  to  2.989. 

Page  34  Paragraph  7.a.iv.(2).  Orange  the  sentence  “For  the  arsenates,  both  CTA  and  KTA 

experience  about  an  8%  loss  while  the  loss  for  RTA  appears  minimal  at  best."  to  “None  of  the 
arsenates  experience  significant  transmission  loss  at  3.17  pm." 

Page  35  Orange  paragraph  7.b.ii.(3)  to  read: 

(3)  3.17/1 .60/1 .064  pm,  III.  While  KTP  has  the  largest  its  absorption  loss  makes 
it  a  poor  candidate.  The  arsenate  KTA  has  the  least  walkoff,  a  slightly  smaller  d^than  KTP  and 
negligible  transmission  loss.  The  only  other  isomorph  candidate,  RTA,  has  only  a  slightly  smaller 
d,,  and  slightly  larger  walkoff  than  KTA.  The  preferred  isomorph  is  KTA. 

Page  37  Oranges  to  Table  12  are  minor  aird  involve  small  changes  to  a  few  values,  all 

Type  I  walkoff  angles,  and  some  changes  in  the  reference  angle  used  to  calculate 
the  relative  walkoff  for  the  other  two  wavelengths.  An  updated  Table  12  has  been 
provided  (note:  k  added  to  updated  table). 

Appendix  B,  Replace  with  the  attached  MathCad  worksheet.  The  original  worksheet  incorrectly 

pages  47-51  treats  the  isomorph  RTA  as  a  negative  crystal  (it  is  positive),  and  incorrectly 
calculates  relative  walkoff.  Additionally,  a  few  other  simplifications  were  made. 
Descriptions  of  the  changes  are  as  follows: 

Eqns  B-3a  and  B-3b  ensure  that  the  crystal  type  is  identified  as  ±1.  B-6a 
and  B-6b  identify  the  correct  angle  0  column  of  the  data  matrix.  Eqn  B-6c 
creates  the  appropriate  phase  match  angle  pair  matrix  for  that  phase  match  type. 
Eqns  B-15a,  B-15b  and  B-15c  correctly  identify  the  a;,  aj,  and  a^,  unit  vectors 
and  eqn  B-16  creates  the  column  vector,  a,!,.  Eqns  B-19a  through  B-19c 
calculate  the  walkoff  angle  for  the  selected  phase  match  type.  Eqns  B-19d  and 
B-19e  setup  for  eqn  B-19g  which  calculates  the  relative  walkoff  angles.  The 
smallest  walkoff  angle  identifies  the  column  which  has  the  correct  relative 
walkoff  angles  (e.g.,  if  pj  is  smaller  than  the  other  two  angles,  then  the  third 
column  has  the  correct  relative  walkoff  angle  values). 

Appendix  D, 

pages  61  -  65 


Attachments: 

1.  Table  12,  Update 

2.  Appendix  B,  Update 

3.  Appendix  D,  Update 


Tables  D-1  through  D-5  are  replaced  by  the  attached  page  of  tables. 


DALE  L.  FENIMORE,  Captain,  USAF 
Electro-Optics  Projects  Engineer 
Electro-Optics  Sources  Branch 


CALCULATE  BiAXIAUUNIAXIAL  CRYSTAL  NONLINEAR  COEFFICIENT  {d,„) 
'  and  WALKOFF  (p)  FOR  POSITIVE/NEGATIVE  CRYSTALS 


OBTAIN  DATA  FROM  DATA  FILES 

M«lrix:=READPRN(pm_d*Urt,)  d  :=READPRN(dijk  «,)  n  :=READPRN(indicesdrt,) 

IDENnFY  ORIGIN  FOR  WORKSHEET  MATRICES 
OUGIN  :=o 


dykrortliecrysUl 


Indices  of  RcflractloD 


0.00 

0.00 

0.00 

0.00 

1.90 

0.00  \ 

0.00 

0.00 

0.00 

3.60 

0.00 

0.00 

,^30 

3.80 

15.80 

0.00 

0.00 

0.00 1  “  "  1 

Select  Positive  or  Negative  Crystal  Type  with  a  positive  or  negative  Number. 


1.7576751  1.7905609 
1.7756109  1.7982079 
1.8278463  1.8655065 


NOTE:  Unless  properly  selecied,  the  calculated  value  for  and  walkolT angles  may  not  be  correct. 
CT:=1  CT:=u{ciS0.1,tf(|CTl>1,CT-lcT'^l.CT))  CT  =  1 


PHASE  MATCH  angles  daunie  LEGEND  [Matrix],  (by  column):  Xi,X2,&)L3;4:  6(1);  6  (0);  8  (111) 
SELECJ  AN  ANGLE  (!)  FOR  WHICH  AN  ANGLE  6  EXISTS  ffor  Tvnr  I  IL  or  IR  Selected  Belowl 


i:=0 


pmid:=2 


X:=[^(Matrix‘^*^^)j  (Matrix*^^  ;  (Matrix^^)ij 


XjI-Xqo  ^2-~Vi 

Xj  =3.176  X  2  =  1-6000 


^3 ‘=^,2 

X3  =  1.0640 


C:=(4  5  6) 


c  :=if(PMIDS3,Co  2.if(PMID=2.Co_^ .Co  o)) 


T  ‘^augment  (Matrix**^ ,  Matrix*^^^) 

Aj  :=(t'^°^)  i'deg  p;  :=(t'’  j  deg 

8:=Ai  (!.;=Pi 

6  =  44.4295*deg  <i  =  0.0000’deg 


SEI JCT  Phase  Match  Category  for  Desired  Nonllnt-ar  Coemdent  fdwrl 
1  -  SHG/SFM  ("Xa")  2  -  DFM-32/OPA-32  ("Xl") 

3  -  DFM-31/OPA-31  ("X2")  4  -  OPO  ("Xl") 

CATEGORY :=4 


(B-la,b,c) 

(8-2) 

1.8026541  \ 

1.8092106 

1.8809392/ 

(B-3a.b) 

(B4) 

(B-5a) 

(B-5b,c,d) 

(B-ia,  b) 

(B-6c) 

(B.7a,b) 
(B-7c.  d) 

(B-8) 


47 


Thb  ngioH  bttint  the  caUuUuionJor  the  tjftdht  noihtinear  eotjpcitia  C4^>- 
Thttt  ralrutalioiu  an  NOT  limited  to  the  prineipU-axis. 


Tbix  is  the  solutioa  for  the  OPTIC  AXIS,  Q 


^  ”2.0  (°1.o)  -(°0.0)  ^,0-18  experiencing. singul«ity. 

["’•“[w'-K/j 


n:=  %L  -Fo.i)  ^,0-18 

P-’L("8.l)'-("o..fJ  J 

r  2  2\A 

^  ("li  -(°o.d 

.  r'-np/-(“o/J  J 


This  uigle,  £i,  is  mevured  from  the  z-wris  to  the  optic  tucis. 
NOTE:  the  10*^8  value  is  intended  to  prevent  unZdErom 


/31.107\ 

Q=  19.167  rdeg 
117.233/ 


nils  Is  used  to  calcuUte  the  poUrizsUon  uigte,  "S".  This  angl*  Is  formsd  by  the  direction  of  the  stow  ray  electric  displacement  unit 
vector,  et,  from  the  surtece  of  the  plane  formed  by  the  wavefront  vector  (E)  and  the  2-axIa 

_ cos(6)-sin(2-^) _ 

cos(e).sin(20) _  /OOOOO  \ 

\2  /„x2  f.2  /,2  Un25  =  0.0000  (B- 

“*(^1.o)  -eos(e)  cosW)  .fonW)  \ 0.0000  / 

cos(o)'sn(2-4i) 

®*dj‘“^(°“(®)‘*“(2<>)®’0.0.'K“»(®)®“(2-4)>0,0,ll))  8^j  =  0.0000  cos(e)sin(2-^)  =0.0000  (B- 

’(““H®o.o)-p5.di)-o.sl 

8:=  (««i(tan26i  (,)^-6^^-0.5  8  =  |  0.000  *deg  (B- 

(«sn(un262^o)  +« «ii)-0-5  '  °°°°  ' 


8  ^j= 0.0000 


(0.000  \ 
0.000  rdeg 
0.000  / 


cos(e)sin(2(|()  =0.0000 


The  following  are  the  solutions  for  the  312)W  (bl)  rays  wd  the  FAST  {b2}  rays.  These  solutions  represent  the  IINTT  VECTORS 
ofthc  electric  dlspUcementvectai  a  Columns  represent  Xy,  Xg,  &X3,  respectively. 

a:=-cos(0).cos((ji)-cos(6(j  Q)-f-$in(^)  sin(SQ_o)  b  :=-co$(e)  sin(4i)  cos^6o_Q) -cos(0)-sin(6Q  q) 

c  I=-cos(6)-cos(ij)-cos  (8y  0)  han('!i)  •sia  (8,  0)  d  :=-cos(e)  an(<^)  cos(8i  0)  - eo$(«^).sin(8y  q) 

e  i  =-cos(6)  •cos(<j»)  'COS  (®2,o)+““W  f: — oos(8}-sin(4i}-cos^82^Q)  — cos(i^)'Sm^62  q) 


(B-na.b) 

(B-nc,d) 

(B-lle,0 


a.  :=^cos(e)  cos((^)-sin^8o  q)  "  *“(4')-'»‘(^,o) 
cc  :=-cos(e).cos(^i).sin(8y  q)  -  sin(^)-cot^6y  q] 
ec  :=-cos(e)-cos(^}>sin^82  q)  - $in(^)-co$^82_Q) 


i(e)'COs(SQ  q]  an(8)-cos(8y  q)  sin(e)-cos^82  q)|  ^i(e)-sin(6Q  q)  sin(e)-sin^8y  q)  sin(e) -501^82  g) 


bb  '• — ^cos(6}-sm(l{i)'sin  (5o,o)  «»(So.o) 

dd  ;~cos(6)-sin(^)-sin  (6i_o)+cos((>)-cos(8i_o) 
Sf  :=-cos(e)  an($)-sin^^  g)  +cos($)-eos(82  g) 


(B-ng.b) 

(B-llij) 

(B-llk.l) 


(B-llm,n) 
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.  ■  Thb  McUon  ciiculttei  the  dettrie  Add's  two  noliriMUon  comiiontnts.  Xi.  X2.  oxl  ^3  («>?  to  bottom),  respccUvely. 

Cohirnn  1  isfcrdieSLOWnyt(0i);  colunm  2  calculates  the  FAST i^s  (<2)- 


('’10,o)^l  (‘*'2,0)^ 

0.5 

(•'20.of,(‘’2i.of  _('>22/ 

0.5" 

.(”0.0)^  (“1.0)^  (“2.0)^. 

.  ("1.0)^  ("2.o)^ . 

[(‘’lo.if ,  (‘•ii.if ,  (’’la.if 

0.5 

,  hl.lf  , 

0.5 

("i.y 

.  (“o.y  ("1.1)^  (“2,1)^ . 

[(‘■o./,  (■'1../,  (“•2/ 

0.5 

(“41./ 1  ‘2 

0.5 

..  (”o,^^  (“2.^*. 

.  (“0,^^  (“1.  )  4  (“2,^^  . 

2 


This  setUon  calnilatK  the  UNIT  VECTORS  for  the  nolarfaaUon  components  for  the  electric  fields,  the  SLOW  rqrs  aic  calculated  by  "al,' 
FAST  rays  1^  '12.'’  The  tows  lepresoit  the  s,  y,  and  z-axis  componenu,  reflectively,  the  columns  rtpitsest  X] .  X2  &  X3  (left  to  right), 
le^ecdvely. 


•>10.0 

*’10.1 

‘’>0.2 

“20.0 

•’20.1 

•’20.2 

^o,o(“o,o)^ 

P^i,o("o,i)^ 

^2,o(°0.^^ 

^o.i(“o,o)^ 

•’’2,l(“0.^^ 

*’11.0 

‘’I1.2 

a2:= 

*^2t.O 

•’21.1 

•’2i,2 

”’0.0' ("1.0)^ 

«>1,o(”l,lf 

^2,o("l,^^ 

"’o.iC^i.o)^ 

'’12.0 

*•*2.1 

**12.2 

•’22.0 

•’22,1 

•’22.2 

.^0,tf(“2,o)^ 

^1,0'(“2,l)^ 

”’2,0' (“2,^^. 

.'^o.iKo)^ 

^1,1-(“2,lf 

*’’2,1' ("2,^^. 

Categories:  CATEGORY  =4  Crystal  Type  tPosilhre  or  Negative 

l=SHG/SFM,2  =  OPO-32/DFM>32  Number  indicates  Type) 

3  =  OPO-31/DFM-31,4  =  OPO 


Identify  which  column  U>  use  for  the  fiecific  phase  match  category  selected 


C  i  :=if(CATEGORY=1,2,if(CATEGORY=2,0,if(CATEGORY=3,1,2))) 
Cj:%CATEGORYai,0,if(CATEGORY=2,1.if(CATEGORY=3,0,0))) 
C  k  :=il(  CATEGORY=1 , 1  ,if(CATEGORY=2,2,  il(  CATEGORY=3,2, 1 ))) 


Setup  the  a,,  aj.  and  aj^  vectors 


Phase  Match  Type 
PMID=2 


Ci=2 

(•i) 

0 

II 

u 

(aj) 

Ck  =  1 

(»k) 

Type  I,  n,  and  in 

a;  :=if(cT>0,a2*^^  '^,al^^  '^) 

a j  :=if(cT>0,if(pMIDS2,a2'^^  ,al''^  ,if{pMlDSS2,al*'^ ,82^^ 

aje  :=if(cT>0,if(pMID=3,a2^^  ,if(pMID=3,nl‘^^ 


(B-12) 


the 


(B-13a,b) 


(B-14a) 

(B-Mb) 

(B-14c) 

(B-15a) 

(B-15b) 

(B-lSc) 
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Setup  the  column-vector 


•jO.o'ko.O 

•jl,0*kl,0 


•jk'l 


•j2.0*''2,0 

•jl.0*'‘2,0+*j2,0*kl.0 


•jO,0*k2.0+*j2,0‘*ko,0 


EffcrtlYt  Wonlln  w  Cwtntlmt 
<>e£f-=l«i-(<*ijk-*jlO| 

WAIXOFF  An(;»l»  (o)  Calculations  (columni ;  SLOW  rays,  column  2:  FAST  rays) 
By  foi7,  starting  at  the  top  ie  the  calculation  for  the  Xi,  Xg.  &  Xj.  respectively 


(B-16) 


(BIT) 


(B-18) 


Sdwt  contcl  yilkcff  ind  reUUvt  wtlkofT  .ceonting  to  ohii»  mitch  category 

p  1  :=Re(if(cT>0.if(PMID=2.Po.vPo.o).K^MID=2,Po,o.Po.l))-v) 
p  2  l=Rc(if(CT^C, if^PMIDSy.Pl  1 , Pi  g]  .if^PMIDS^.Pl  Q.Pl , l))"'') 
p  3  :=Re(if(cT>0.p2j.P2,o)->'} 

/I  2  2\ 

T:=^  1  2|  t:=if(pMIDS3,T'^^.if(pMID=2,T‘^^^,f'°^)) 


(B-19a} 

(B-19b) 

(B-19c) 

(B-19d) 


(B-19e,  f) 


P2 


n 

0.18 


P3 


Create  the  Matris  for  Writing  a  Une  of  Data  to  the  Data  Rl> 
Writes  Data  to  the  Data  RIe  Specified 


<tau:=(Xi  Xj  X3  e-v  (It-v  deff  Pi  P2  P3) 
APPENDPRN(dw_info  :=(UtM 


(B-19g) 


(B-2O1) 

(B-20b) 
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Table  D-1.  Nonlinear  Characteristics,  KTP 


XtAj/X] 

Type 

t  m) 

PtWPa 

1.064/1.064/0.532 

1 

45.64,36.93 

47.63/47.63/*8.49 

0.423 

li 

90.00,25.99 

3.47/*0.0/4.87 

3.268 

1.35/1.35rt).675 

i 

3724,31.27 

44.53/44/53^8.16 

0.366 

II 

58.70,0.00 

*0.0/44.68^.0 

2.989 

4.043/1.444/1.064 

1 

36.31,31.61 

38.07/43.92^7.12 

0.438 

II 

44.33, 0.00 

*0.0/48.86/0.0 

2.387 

3.176/1.60/1.064 

1 

34.81,30.10 

39.62/42.70/*7.32 

0.398 

II 

44.85, 0.00 

*0.0/48.55«).0 

2.412 

III 

77.36,0.00 

20.13/*0.0/0.0 

3.496 

Table  D-3.  Nonlinear  Characteristics,  KTA 


1.064/1.064/0.532 

1 

48.64, 39.89 

45.52/45.52/^.26 

0231 

1.35/1.35A).675 

1 

38.62, 34.09 

43.13/43.13^5.31 

0.164 

II 

61.67,0.00 

*0.0/39.43/0.0 

2.781 

4.043/1.444/1.064 

1 

33.31,64.32 

31.41/39.0^*3.96 

0.406 

II 

45.72, 0.00 

0.0/45.68^0.0 

2207 

3.176/1.60/1.064 

1 

32.41,  53.98 

34.73fl827/*5.08 

0.174 

11 

46.01,0.00 

*0.0/45.10/0.0 

2.220 

III 

90.00, 16.40 

*0.0/1.19/1.66 

3.154 

Table  0-5.  Nonlinear  Characteristtcs,  CTA 


1.064/1.064/0.532 

1 

64.20, 41.82 

24.4524.58^9.65 

0.217 

1.35/1.35/0.675 

1 

II 

51.56,27.57 
90.00, 53.40 

31.6321.63/*9.06 

7.80ro.0/7.60 

0.226 

2.492 

4.043/1.444/1.064 

1 

II 

34.16, 61.88 
61.51,0.00 

21.2724.47/*8.88 

0.0/30.83/*0.0 

0.135 

3.110 

3.176/1.60/1.064 

1 

1! 

31.97, 59.41 
62.60, 0.00 

20.9822.79r9.66 

•0.029.85/0.0 

0.177 

3.145 

Table  D-2.  Nonlinear  Characteristics,  RTP 


(0.1>) 

PiWp3 

1.064/1 .06441.532 

52.08,23.32 
90.00, 41.70 

48.48/48.48ri3.62 

11.55r0.0/13.09 

0.128 

2.827 

1.35/1.35rt).675 

26.89, 53.47 
65.95, 0.00 

27.1527.15ri7.37 

0.0/40.83r0.0 

0.324 

3.225 

4.043/1.444h.064 

27.81, 5238 
5279, 0.00 

21.3528.33ri6.87 

*0.051.38/0.0 

0.355 

2.751 

3.176/1.60/1.064 

25.94,49.60 
53.16, 0.00 
90.00,45.35 

20.9824.62ri7.99 

*0.0/50.74/0.0 

*0.0/11.33/11.60 

0.416 

2.768 

2.720 

Table  IM.  Nonlinear  Characteristics,  RTA 


1.064/1.0645.532 

53.11,41.37 

38.9228.99r6.98 

0.469 

1.35/1.355.675 

41.57,36.72 

38.7228.72r6.60 

0.496 

70.40.0.00 

*0.027.10/0.0 

3.357 

4.043/1.444/1.064 

40.48, 13.55 

26.22/40.05r2.89 

1.465 

4248, 0.OC 

0.0/40.95r0.0 

2320 

3.176/1.60/1.064 

37.36, 25.56 

31.0927.54r5.15 

0.952 

44.43, 0.00 

0.0/40.92r0.0 

2.413 

90.00. 58.63 

*0.02.782.22 

2355 

Notes; 

•  Wavelengths  given  in  mid 

•  k'Vector  {^ven  in  degrees 

•  Walkoff  (p)  given  In  mllllradans.  Starred  angle  is  the 
reference  angle  from  which  the  other  two  angles  are 
measured 

•  d^  given  In  terms  of  pm/V 


These  tables  replace  those  found  on  pages  61  -  65  of  the  final  report. 
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DEPARTMENT  OF  THE  AIR  FORCE 

WRIGHT  LABORATORY  (AFMC) 
WRIGHT-PATTERSON  AIR  FORCE  BASE.  OHIO 


21  March  1994 


MEMORANDUM  FOR  WL/DOA  BLD  22 

\jL/prx^ 


FROM:  WL/ELO  Bldg  22B 
2700  D  St  Ste  2 

Wright-Patterson  AFB  OH  45433-7405 

^'02.9. 

SUBJECT:  ChangesAlpdates  to  WL-TR-9355093,  Three-Wave  Nonlinear  Interactions  in  KTi0P04  (KTP) 
and  KTO  Isomorphs 


As  with  many  publications,  there  are  inevitable  errors  in  the  final  product  which  need  correction— and  this 
report  is  no  exception.  Following  are  changes  and  corrections: 

Note  The  distinction  between  “r^”  and  “wave”  is  somewhat  loosely  treated  in  the  report. 

Propeily  considered,  treat  k  ,3  as  a  reference  to  a  “wave”  and  as  a  reference 
to  a  “ray.”  Specifically,  refers  to  the.  wave  nornial  vector  (direction  of 
propagation)  and  S^  refers  to  the  ray  vector  (direction  of  energy  flow). 

Page  4  Figure  3  incorrectly  shows  the  angle  <j>  as  being  referenced  from  the  y-axis.  Modify 

the  figure  to  indicate  reference  as  being  fiom  the  x-axis. 

Page  6  Figure  4  was  obtained  from  reference  [4].  Credit  for  this  figure  was  inadvertently 

removed  from  the  Final  Report. 

Page  10  Last  paragraph.  Replace  the  sentence 

The  value  of  d,„  is  dependent  upon  the  direction  of  the  t...  we  will  use  the  electric 
displacement  vectors,  to  c^culate  a,,  a,,  and  a*). 

with 

The  value  of  dependent  upon  We  do,  however,  knowthat’Sli?  and  that^l^.  We 

also  know  thatD  and^  are  related  in  terms  of  the  polarizability  tensor.  Therefore,  we  can 
use  the  electric  displacement  vectors, '5, (coj,  to  calculate  dj,  and  a^  of  eqn  16. 

Page  12  Change  the  indices  of  refraction  subscripts  in  eqn  23  from  “1,2  and  2”  to  “x,  y,  and 

z”,  respectively. 

Page  14  2nd  Line.  Change  “...the  same  as  {^ijcoj}  (©„))."  to  “...the  same  as  {a^2(<^J} 

Paragraph  3.b.ii,  2nd  line.  Change  “...which  exactly..."  to  “...which  are  exactly..." 

Page  15  Paragraph  3.b.vii.(l).  Change  “...exactly  the  value"  to  “...exactly  the  same  value" 


Page  19 


Table  2,  last  cell.  The  Sellmeier  equation  should  read: 

n?  =  Ai  +  (Bi/[l-(CA)"]-DiX^ 

Page  30  Table  11.  Change  d^f  for  KTP,  0.675  pm-SHG,  Type  II  from  2.603  to  2.989. 

Page  34  Paragraph  7.a.iv.(2).  Change  the  sentence  “For  the  arsenates,  both  CTA  and  KTA 

experience  about  an  8%  loss  while  the  loss  for  RTA  appears  minimal  at  best."  to  “None  of  the 
arsenates  experience  significant  transmission  loss  at  3.17  pm." 

Page  35  Change  paragraph  7.b.ii.(3)  to  read: 

(3)  3.17/1 .60/1 .064  pm,  Type  III.  While  KTP  has  the  largest  d^,  its  absorption  loss  makes 
it  a  poor  candidate.  The  arsenate  KTA  has  the  least  walkoff,  a  slightly  smaller  d„than  KTP  and 
negligible  transmission  loss.  The  only  other  isomorph  candidate,  RTA,  has  only  a  slightly  smaller 
da  and  slightly  larger  walkoff  than  KTA.  The  preferred  isomorph  is  KT/t. 

Page  37  Changes  to  Table  12  are  minor  and  involve  small  changes  to  a  few  dj,,  values,  all 

Type  I  walkoff  angles,  and  some  changes  in  the  reference  angle  used  to  calculate 
the  relative  walkoff  for  the  other  two  wavelengths.  Aji  updated  Table  12  has  been 
provided  (note:  k  added  to  updated  table). 

Replace  with  the  attached  MathCad  woiksheeL  The  original  worksheet  incorrectly 
treats  the  isomorph  RTA  as  a  negative  crystal  (it  is  positive),  and  incorrectly 
calculates  relative  walkoff.  Additionally,  a  few  other  simplifications  were  made. 
Descriptions  of  the  changes  are  as  follows: 

Eqns  B-3a  and  B-3b  ensure  that  the  crystal  type  is  identified  as  ±1.  Etpis  B-6a 
and  B*6b  identify  the  correct  angle  0  colutrm  of  the  data  matrix.  B-6c 
creates  the  appropriate  phase  match  angle  pair  matrix  for  that  phase  match  type. 
Eqns  B-15a,  B-15b  and  B-15c  correctly  identify  the  aj,  aj,  and  a^  unit  vectors 
and  eqn  B-16  creates  the  column  vector,  aj^.  Eqns  B-19a  through  B-19c 
calculate  the  walkoff  angle  for  the  selected  phase  match  type.  Etps  B-19d  and 
B-19e  setup  for  eqn  B-19g  which  calculates  the  relative  walkoff  angles.  The 
smallest  walkoff  angle  identifies  the  column  which  has  the  correct  relative 
walkoff  angles  (e.g.,  if  p3  is  smaller  than  the  other  two  angles,  then  the  third 
colurrm  has  the  correct  relative  walkoff  angle  values). 

Appendix  D,  Tables  D-1  through  D-5  are  .^placed  by  the  attached  page  of  tables, 
pages  61  -  65 


DALE  L.  FENIMORE,  Captain,  USAF 
Electro-Optics  Projects  Engineer 
Electro-Optics  Sources  Branch 

Attachments: 

1.  Table  12,  Update 

2.  Ai^ndix  B,  Update 

3.  Ajqrendix  D,  Update 


Appendix  B, 
pages  47  -  51 


Table  D-1.  Nonlinear  Characteristics,  KTP 


Table  D-2.  Nonlinear  Characteristics,  RTP 


1.064/1.064W.532 


1.35/1 .35«).675 


4.043/1.444rt.064 


3.178/1.60/1.064 


I 


PilpJPi 


45.64,36.93  47.63/47.63/*8.49 
90.00,25.99  3.47/*0.0/4.87 


37.24,3127  44.53/44^3^8.16 
58.70, 0.00  *0.0/44.68/0.0 


36.31, 31.61  38.07/43.92/7.12 
44.33, 0.00  *0.0/48.86A).0 


34.81,30.10  39.62/42.70/7.32 
44.85, 0.00  *0.0/48.55ra.0 
77.36, 0.00  20.13/*0.0/0.0 


Table  IM.  Nonlinear  Characteristics,  KTA 


1.064/1.0642.532 

m 

48.64, 39.89 

45.52/4522^526 

0231 

1.35/125/0.675 

■ 

38.62,34.09 

43.13/43.13^5.31 

0.164 

61.67, 0.00 

*0.029.432.0 

2.781 

4.043/1.444^.064 

1 

33.31, 64.32 

31.4129.05^3.96 

0.406 

II 

45.72,0.00 

0.0/45.68r0.0 

2207 

3.176/1.60/1.064 

JM 

32.41, 53.98 

34.732827/*5.08 

0.174 

46.01,0.00 

*0.0/45.102.0 

2220 

90.00, 16.40 

*0.0/1.19/1.66 

3.154 

1.064/1.064 

12.532 

1.35/1.352 

.675 

4.043/1.444 

1/1.064 

3.176/1.60/ 

1.064 

Table  D-5.  Nonlinear  Characteristics,  CTA 


6420,41.82  24.4524.58^9.65 


51.56,27.57  31.6321.63^9.06 
90.00,53.40  7.80^0.0/7.60 


34.16, 61.88  212724.47/*8.88 
61.51,0.00  0.0/30.83^0.0 


31.97,59.41  20.9822.79/*9.66 
6^60, 0.00  *0.0/29.85/0.0 


X,/Xj/X3 

PJPJPZ 

1.064/1.0642.532 

52.08,23.32 
90.00, 41.70 

48.48/48.48/*13.62 

11.55/*0.0/13.09 

1.35/1.352.675 

26.89, 53.47 
65.95, 0.00 

27.1527.15/*17.37 

0.0/40.83^0.0 

4.043/1.444/1.064 

27.81, 5238 
5279, 0.00 

21.35/28.33n6.87 

*0.021.38/0.0 

3.176/1.60/1.064 

25.94, 49.60 
53.16,0.00 
90.00, 45.35 

20.9824.62/*17.99 

*0.020.742.0 

*0.0/11.33/11.60 

Table  0-4.  Nonlinear  Characteristics,  RTA 


1.064/1.0642.532 

53.11,41.37 

38.92/38.99/*6.98 

1.35/1.352.675 

41.57,36.72 

70.40.0.00 

38.7228.72/*6.60 

*0.0/27.102.0 

4.043/1.444/1.064 

40.48, 13.55 
4248, 0.00 

26.22/40.05^289 

0.0/40.95/*0.0 

3.176/1.60/1.064 

37.36,25.56 
44.43, 0.00 
90.00, 58.63 

31.0927.54/*5.15 

0.0/40.92/*0.0 

*0.02.782.22 

Notes: 

•  Wavelengths  (X„X2,X,}  given  In  pm 

•  k'Vector  given  in  degrees 

•  Walkofi  (p)  given  in  millitadians.  Starred  angle  is  the 
reference  angle  from  which  the  other  two  angles  are 
measured 

•  d^  given  in  terms  of  pm/V 


These  tables  replace  those  found  on  pages  61  -  65  of  the  final  report 


CALCULATE  BIAXIAUUNIAXIAL  CRYSTAL  NONLINEAR  COEFRCIENT  (</,«) 
and  WALKOFF(p)  FOR  POSmVE/NEGATlVE  CRYSTALS  _ 

OBTAIN  DATA  FWOM  DATA  Fn.ES 

MtirU  :=READPRN(pm_<i»i  «•)  ijk  :=READPRN(<lijk  «,)  n  :=READPRN(indicesd 

IDENTIFY  ORIGIN  FOR  WORKSHEET  MATRICES 
ORIGIN  :=o 

dykforthecrjrsUI  Indices  of  Rdt^Uon 


10.00 

0.00 

0.00 

0.00 

1.90 

1.7905609 

1  / 1.7576751 

dijk=  0.00 

0.00 

0.00 

3.60 

0.00 

0.00  1 

I  n=  1.7756109 

1.7982079 

\2.30 

3.80 

16.80 

0.00 

0.00 

0-00  /  1  ^ 

1.8655065 

\  1.8278463 

Sflect  Poritivg  nr  Negative  Crystal  Type  wiih«  POSITIVE  or  NEGATIVE  Number. 

NOTE:  Unless  properly  selected,  the  odculaled  value  f°r<f^  and  walkolT  angles  may  not  be  correct. 

cr:=i  cT;=if(cpso.i  if(lcrl>i.cT  |cr'1.cT))  cr=i 

PHASE  MATCH  angles  date  die  LEGEND  (Matrix],  (by  column):  Xi,X2,&X3;^:  6(0: 6  (B);  6  (DO 
SELECT  AN  ANGLE  «1)  FOR  WHICH  AN  ANGLE  6  EXISTS  (for  Type  1.  H.  or  IE  Selected  Belowl 


i:=o 


pmid:=2 


X:=[(Matr«^**^)j  (Matrix*^^^)}  (Matrix*^^);] 

^2 -“Vi  ^3‘“V2 

Xi=3.176  X  2  =  1.6000  X  3  =  1.0640 

C:=(4  5  6)  c:=if(PMII>=a.Co  2.i((PM®=2.Cg 


T  I-wgmeniVMatrix^^ 


.Matrix*^^^) 

Pi:=(T‘=’^)ideg 


o:=Ai 


H-Pi 


e=44.4235*deg  ^  =  0.0000*deg 


?5FT.ErT  Wiaee  Match  Category  for  Desired  Nonlinear  Codndtnt  fdeiTl 
1  -  SHG/SFM  r\3”)  2  -  DFM-32/OPA-32  ("Xl") 

3  -  DFM-31/OPA-31  ("X2")  4  -  OPO  ("Xl") 


(B-la,  b,  c) 

(B-2) 

1.8026541  \ 

1.8092106 

1.8809392/ 

(B-3a,b) 

(B4) 

(B-5a) 

(B-5b,c,d) 

(B-6a,b) 

(B-6c) 

(B-7a,b) 

{B-7c,d) 


category:  =4 


(B-8) 


This  ntiom  b*giiu  tkt  caleuUulon/or  Iht  ejfteliv*  noH-Unear  eoeXfitUat  (i^)- 
Thttt  eubutallons  art  NOT  UmUtd  to  the  prineipU-axis. 


This  is  the  solution  for  the  OPTIC  AXIS,  Q 
2  2*  *5 

("i.rf  -Ka) 

2  2" 

Q:=  «ia  -  ^.10-18 

[“l.l[(„2.l)"-Kin 

2  2* 

,a,  !M.  (°ti  -(*oi  ^,0-18 

,  h-nw^-Kil  J 


This  angle,  O,  is  measuicd  firom  the  s-axis  to  the  optic  axis. 
NOTE:  die  10  value  is  intended  to  prevent  tanZS  firom 
experiencing  a  singularity. 


'31.107  \ 

Q=  19.167  ‘deg 
,17.333/ 


TUt  is  used  to  calculate  the  polarization  angle,  "S”.  This  sngla  is  formed  by  the  diractlon  of  ths  slow  ray  sisctric  displacoment  unit 
vector,  e^,  from  the  surtax  of  the  plane  formed  by  the  wavefront  vector  (A)  and  the  z-exis 

r _ coi(e)-tin(2-4i) _ 1 

\2  .  /„\2  /^\2  /.\2  .  /,x2 


tan2Sl= 


_ cos(6)-sin(2-0) _ 

cos(6)-an(2'^) 

9  9  9  9  C 

cot^Q2,o)  •sin(6)  —  cos(6)  •cos(^)  ■t-sin(4>) 


/  0.0000 
tan2S=  I  0.0000 


^j;=if(ooc(e)-sin(2-^)aO.O,if(cos(0)-sin(2-^)>O,O.Jt)) 

8  ^j  =  0.0000 

(alan(tan2^_g)4.S^^.0.5' 

1 0.000  \ 

:=  (alaa(tan2a,  g)  4-8 

8  =  1  0.000  ♦« 

(alan(tan282  g)-i-8^^.0.5 

\  0.000/ 

cos  (e)-sin(24)  =0.0000 


The  following  arc  the  solutions  for  the  SLOW  (bl)  rays  and  the  FAST  (b2)raj^  These  soliitionsrqitcsaittbellNITVECTORS 
of  the  dectric  tflsplacement  vectors.  Columns  represeat  Xy,  Xg,  8c  Xy,  respectively. 

a:=-cos(6)-cos(4)-co$^6g^(lj)  .f-sin(^)<sin^SQ^Q)  b  :=-cos(8)'Sin(^)'Cos  (®o,o)  -cos(^)-  nn(5o,o) 

c  •^cos(6)*cos(^)<cx)s  (Sf  Q)+an((!i)-sin(6i  0)  d:~cos(e)  siD((!i)-eof(6i  q)-eos(4>).sin(8i  q) 

c ;  ^  coi(0)  <cx>s(^)  'COS  (82,0) +*“(♦)  *&a  (^.0)  r:— CO*  (e)  •sia(4i)-  cos(82  0)  -cos(4i)-sin(82  0) 


^cos(e)-cos(^)-sin(8Q  g)  - sin(4))-cos^So,Q) 

=^cas(e)-cos(^-sia^8.|  g)  -  sin(4i)-cos^8^  g) 
=-cos(e)-cas(^)-sin^82  g)  ~  ”“(^)’““(*i2,o) 


bb : — ^oos(o)'sin(^)'sn  (80,0)  +«»»(<')-=«(Sg_o) 
dd  l^cos(6)-sin((^}-sui  (81  g)+cos(4>)-cos(8i_o) 
fif  :=-c«(e)  •sin(^)  -sin^^^  g)  -f-cosC^}  -cos  (82^  g) 


^iin(6)'Cos^8g^gj  sin(6)-cos^8^^g^  an(o)  •005(82  0) 


pn(e)-sin(5g  g)  sin(e)-sin(8,,  g)  sin(e}-sin(S2  g) 


(B-lla,b) 

(B-llc,d) 

{B-lle,0 


(B-llgJi) 

(B-llij) 

(B-llk,!) 


(B>llm,n) 
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Thh  section  r»lcul»tM  the  dwtric  Held's  two  poIitImUoii  comnontnls.  X] ,  X2.  and  X3  (top  to  bottom),  respectively. 
Column  1  is  for  the  SLOW  rays  (^1);  column  2  calculates  the  FAST  rays  (^2). 


< 


(‘’10,o)^  1  {*’»1,o)^  1  (*’l2,o) 

0.5 

0.5‘ 

.  (“0,o)^  (”l,o)^  (°2,o) 

('’l0,l)^  1  (‘’ll.l)  1  ('’l2,l) 

0.5 

(“O.o)^  ("l.o)*  (“2.o)^ , 

(‘20,1)"  1  {"21.1)",  ('’I2.1)" 

0.5 

.  (“O,!)"*  (“l.l)"*  (“Z.l)"*. 

0.5 

.(“0,/  ("i.y 
('’2o.y|  '*2  2^^  ('’22.2)^ 

0.5 

1  (”0,^^  (“1,^^  (“2.^'*. 

.(“0,2)"  (“1,)  4(“2,2)^ 

2 


This  sedlon  calfulatw  th>  TINIT  VECTORS  for  the  polarization  components  for  the  electric  Helds.  The  SLOW  rays  are  calculated  by  "al the 
FAST  ri^s  1^  "a2.''  The  rows  iqxescat  the  s,  >,  and  z*axis  components,  reflectively,  the  columns  represent  X] ,  X2  &  X3  (left  to  right), 
respectively. 


*’10,0 

‘’I0.I 

‘’10.2 

2 

*’’o,o(”o,o) 

^2,0‘(“0,^^ 

•’M.o 

‘’11,1 

•’I1.2 

*^o.tf("i,o)^ 

«’i.o(»i,y 

‘’’2,o'(“i,y 

‘’»2.0 

'’12.1 

'’12.2 

>.0-(»2,0)" 

”’i,o(“2.y 

‘*2,0'(“2.^^. 

'’20.0 

‘’20.1 

”20.2 

‘’•’o.i'N.o)^ 

‘^l.l-Klf 

‘’’2,1‘(”0,^^ 

'’2l.0 

‘’21.1 

”21.2 

‘^o.i-(“i.of 

«’l.1-("l.lf 

‘’”2.i("i.y 

‘’22.0 

”22.1 

”22.2 

.^0.1'(“2.o)^ 

‘^l.l-Klf 

‘’’2,1' ("2.^^. 

CatgaricF  category = 4 

1  =  SHG/SFM,  2 = OPO-32/DFI4-32 
3  =  OPO-31/DFM-31, 4  =  OPO 


rrvstalTvne  (Positive  or  Negative 
Number  indicates  Type) 

CT  =  1 


Identify  which  column  to  use  for  the  fiecific  phase  match  categoiy  selected 
C  i  :=if(CATEGORY=1,2,if(CATEGORYSS2,0,if(CATEGORY=3.1,2))) 
C  j  :=if(CATEGORY=1 ,0,if(CATEGORiTO,1.if(CATEGORY=3,0,0))) 
C  k  :=if(  CATEGORY=1 , 1 ,  if(  CATEGORY=2, 2,  if(  CATEGORY=3, 2,1))) 


Phase  Match  Type 
PMID=2 


Ci=2 

(8i) 

(B-14a) 

Cj=o 

(aj) 

(B-14b) 

Ck  =  1 

(•k) 

(B-14c) 

TVpeLILnndni 

a  j  :=if(cT>0,a2‘^^  ,al^^  ‘^) 

a  j  :=:f  (cT>0,  ifG‘MIDS2,  a2''‘^  ,  al"^^  ,  if(pMID=2,  al^^  ,  a2’'^ 

uk  :=if(cr>0,if(pMID=3,a2^*^  .al"'^  ,if(pMID=3,al'^^  ''^)) 


(B-15a) 

(B-15b) 

(B-lSc) 
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CONSOLIDATED  INFORMATION  Crystal  rta 


Electric  Displacement  Unit  Vectors,  bl  =  slow  b2  =  fist  r«y 


Polmzation  Unit  vectors,  al  =  "slow"  rays,  a2  =  "fast"  rqrs 


^.7409  -0.7421 

-0.7432^ 

,  /  0.0000 

0.0000 

0.0000  ' 

al  s| 

0.0000  0.0000 

0.0000 

a2°  1.0000 

1.0000 

1.0000 

i  0.6716  0.6702 

0.6691  / 

'  \  0.0000 

0.0000 

0.0000  i 

O.OCOO  0.0000  0.0000 
1.0000  1.0000  1.0000 
0.0000  0.0000  0.0000 


G.7141  -0.7141  -0.7V' l\ 
[)000  0.0000  0.0000 
7000  0.7000  0.7000/ 


Polarization  Unit  vectors 

aj  =  "Generated  Wavelength",  aj  and  a|c  =  "Applied"  Wavelengths 


/  0.0000  \ 

/  0.0000  \ 

^0.7421  \ 

ai  =  ’  ".0000 

aj=  1.0000 

aj^s:  I  0.0000  1 

\  00000/ 

\  0.0000/ 

\  0.6702/ 

ajaic  Cohir :  i  vector  => 


Polarization  Components  for  the 
Electric  Fields  (col  =  slow  &  fast 


■  0.0000  ■ 

respectively) 

0.0000 

0.0000 

/  0.3120  0.3172^ 

0.6702 

Pb  =  0.3001  0.3093 

0.0000 

\oE957  0.3055/ 

.-0.7421  J 

Miscellaneous  Information 

dlJTer  . 

i 

/o.;o 

0.00 

V  ■> 

0 

0 

c 

1.90 

It 

•0 

000 

0.00 

0.00 

3.60 

0.00 

i2.30 

3.80 

15.80 

0.00 

0.00 

Indices  of  Reft-adlon  (X  by  columns) 

0.00  \  / 1.7577  1.7906  1.8027  \ 

0.00  n=  1.7756  1.7982  1.8092 
0.00/  \l.f:78  1.8655  1.8809/ 


Polarization  Veciu. 
d_ajak  :=d  ijf*  jk 


WalltoB’,  column  I  =  slow  rays 
column  2 = fast  rays,  X  by  rows) 
in  degrees  and  milliradians,  resp. 


f  0.0000^ 

, 

[22386 

12074*10”®i 

d_ajak  =  | 

2.4128 

i  0.0000  j 

1 

1  R3443 
iR4302 

12074*10r®i 

0.0000 

Nonlinet.*  Coefficient  (pmiv) 

deff  =  R4128 


p-10^  = 


39.0683 

40.9163 

42.4158 


2.1073*10 

2.1073*10~®i 

0.0000 


Relative  Walkoif  angle  in  mUUradians 

Use  associated  walkoff  column  below  for  the  smallest  walkoff  angle 


8  =  44.4295 ‘deg  ^  =  0.0000‘deg 

X  1=3.1761  X2  =  1.6000 

X3  =  1.0640 

CATEGORY  =4 
1-  SHG/SFM,  20PA/DFM-32 
3.0PA/DFM-31,4-0P0 

Crystal  Type 

Phase  Match  (4/.„umber 

indicates  Qpe) 

PMiD=2  cr  =  i 

Angleof  /31.1070  \ 

OPnCAxis  Q=  19.1667  'deg 
\  17.3332/ 


Polarization 

angle 


(0.0000  \ 
0.0000  ‘deg 
0.0000/ 


(Pl>  P2i  tv  P3)  at  indicated  to  the  tight. 


PI 

P2 

P3 

0 

40.9163 

0 

40.9163 

40.9163 

40.9163 

0 

40.9163 

0 

Worksheet  Walkoff  values  (by  selected  phase 
match  Qpe)  in  degiees 

Pi=o 

P2  =  R344 


P3=0 
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DEPARTMENT  OF  THE  AIR  FORCE 

WRIGHT  LABORATORY  (AFMC) 
WRIGHT-PATTERSON  AIR  FORCE  BASE.  OHIO 


FROM: 

SUBJ: 

TO: 


WL/DOA 

Wright-Patterson  AFB  OH  45A33-6523 


Notice  of  Changes  in  Technical  Report (s) 


Defense  Technical  Information  Center 
Attn:  DTIC-OCC 
Cameron  Station 
Alexandria  VA  22304-6145 


Please  change  subject  report (s)  as  follows: 


JUDY  K.  SZCZUR 

Chief,  STINFO  and  Technical  Editing  Divison 
Operations  and  Support  Directorate 


